BBS  file  .corf 


r~, 

'v»\- 


Real-Time  and  Applications  Software  for  the 
Neil  Brown/W.H.O.I.  CTD  Microprofiler 


Lawrence  J.  Rose.nbl.um 
Oceanographic  Computer  Applications  Group 


December  30,  1980 


NAVAL  RESEARCH  LABORATORY 
Washington,  D.C. 


Approved  lor  pi  ilie  release  disirihuiion  unlimileil 


security  classification  of  this  page.  Data  tnt«r*dj 


REPORT  DOCUMENTATION  PAGE 


RE  PORT  NUMBER 


NRL  Report  8444 


4  TiTlE  (and  Subtitle) 


l2  GOVT  ACCESSION  NO.  3 

Sd  ~4c  ''V  <■> 


jtUAL-JIMK  AND  APPLICATIONS  JSOFTWARL  FOR  THE 

NFIL  BROWN/W.H.O.I.  CTD  MICROPROFI LLR  . 

-  * 


READ  INSTRUCTIONS 
BEFORE  COMPI.KTINO  FORM 


RECIPIENT'S  CATALOG  NUMBER 


TYPE  OF  REPORT  A  PERIOD  COVERED 


[Final  report 


PERFORMING  ORG  REPORT  NUMBER 


1  Lawrence  J./Ros 


Rosenblum 


CONTRACT  OR  GRANT  NUM0ERr«J 


9  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


Naval  Research  Laboratory 
Washington,  DC  20375 


program  Element  project,  task 
AREA  a  WORK  unit  numbers 


J.  O.  68056-6013 


11  CONTROLLING  OFFICE  NAME  AND  ADDRESS  12  REPORT  DATE 

j  :  December  30,  1980 

1  ~  NUMBER  OF  PAGES 

46 


14  MONITORING  AGENCY  NAME  a  ADDRESSf/f  dlttaram  from  Controlling  Oftlca )  15  SECURITY  CLASS  (ot  thta  raporl) 

UNCLASSIFIED 


DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16  DISTRIBUTION  STATEMENT  (ol  this  Report) 


Approved  for  public  release;  distribution  unlimited. 


17  DISTRIBUTION  ST  ATEMENT  (ol  the  abatracl  antarad  In  Block  20,  II  dlllarant  Irom  Raport) 


19  KEY  WORDS  f Contlnua  on  ravaraa  alda  If  nacaaamry  and  Idantlly  by  block  numbar) 


Oceanography 

Ocean  microstructure  measurements 
CTD  microprofiler  data  acquisition 
CTD  microprofiler  real-time  data  display 


Real-time  processing 
Minicomputer  applications 


20  ABSTRACT  f  Continue  on  ravaraa  alda  H  nacaaamry  and  Idantlly  by  block  numbar) 


~->’The  Neil  Brown/W.ll.O.I.  CTD  Microprofiler  has  been  extensively  used  by  a  variety  of  NRL 
researchers  in  at-sea  experiments.  This  paper  reports  on  software  for  both  real-time  data  acquisition 
and  display  and  for  processing  collected  data. 


CONTENTS 


THE  CTO  MICROPROFILER . 

Introduction . 

Interfacing . 

REAL-TIME  COLLECTION  AND  DISPLAY  SOFTWARE 

The  Real-Time  Program  Set . 

System  Requirements . 

Data  File  Structure . 

Data  Acquisition  and  Display . 

Real-Time  Programs— FORTRAN  Source  Listings . 

USER’S  GUIDE  TO  REAL-TIME  CTD  PROGRAMS . 

Preparation . 

Program  Start-Up . 

Special  Functions . 

DATA  PROCESSING  PROGRAMS  . 

Parameters  of  Interest . 

User’s  Guide  to  Program  PLTBK . 

Examples  of  Program  PLTBK  Graphics . 

Data  Processing  Programs  — FORTRAN  Source  Listings 
REFERENCES . 


APPENDIX— Scientific  Subroutines— FORTRAN  Source  Listi 


1 

I 

I 

3 

3 

3 

4 

5 
9 

24 

24 

24 

26 

27 

27 

28 
28 
32 

39 

40 


REAL-TIME  AND  APPLICATIONS  SOFTWARE  FOR  THE 
NEIL  BROWN/W.H.O.I.  CTD  MICROPROFILER 

THE  CTD  MICROPROFILER 

Introduction 

The  W. HOI. /Brown  conductivity,  temperature,  and  depth  (CTD)  profiler  is  a  device  to  make 
fine-scale  measurements  of  these  physical  parameters  in  the  ocean.  The  CTD  was  first  used  by  the 
Naval  Research  Laboratory  (NRL)  in  July  1977  and  has  been  extensively  used  in  field  work  in  physical 
oceanography,  chemical  oceanography,  acoustics,  and  special  projects  ever  since.  In  August  1979  one 
of  the  NRL  CTD  devices  was  modified  to  include  a  fluorometer.  All  field  work  utilizing  the  CTDs  has 
used  on-board  computer  systems  of  the  Oceanographic  Computer  Applications  Group.  This  report 
documents  software  written  by  the  author  for  the  real-time  acquisition,  storage,  and  display  of  CTD 
data  as  well  as  application  programs  for  processing  collected  data. 

The  CTD  consists  of  a  shipboard  data  terminal  deck  unit  and  an  underwatei  unit  which  provides 
continuous  sampling  of  the  variables  as  it  is  lowered  through  the  water.  Conductivity  is  measured  with 
a  miniature  four-electrode  ceramic  cell.  Temperature  is  sensed  by  electronically  combining  the  outputs 
of  a  high-speed  thermistor  and  a  platinum  resistance  thermometer.  Pressure  is  sensed  by  a  strain-gauge 
pressure  transducer.  Data  from  the  underwater  unit  is  transmitted  in  real  time  to  the  shipboard  data 
terminal  through  a  single  electrical  conductor  in  the  core  of  the  steel  cable  used  to  support  the  instru¬ 
ment  in  the  water.  The  data  are  in  Teletype  format  using  frequency-shift  key  modulation  of  5-  and  10- 
kHz  uc  signals  superimposed  on  the  dc  power  supply  to  the  underwater  unit  through  the  same  cable. 
The  deck  unit  decodes  the  signals,  provides  digital  data  in  both  parallel  and  serial  format  for  computer 
processing,  and  displays  conductivity,  temperature,  and  pressure.  For  greater  detail  on  the  CTD  see 
references  1  to  V 

The  NRL  Oceanographic  Computer  Applications  Group  uses  Hewlett-Packard  (HP)  1000-series 
computers  and  a  variety  of  peripherals  depending  on  individual  requirements.  For  the 
September/October  1979  cruise  aboard  the  USNS  Lynch  two  HP21MX-E  computers  were  used,  one  for 
taking  real-time  data  from  various  physical  oceanographic  data  collection  devices  and  ship's  sensors  and 
the  other  for  recording  CTD  and  XBT  data.  The  CTD  system  produced  real-time  plots  on  a  HP2647 
terminal  and  used  an  HP2631G  line  printer  for  producing  hard  copy  of  these  plots.  Minor  program 
modifications  are  necessary  for  other  peripherals.  Figure  1,  taken  aboard  the  USNS  Hayes  in  1977, 
shows  a  typical  shipboard  computer  system. 

Interfacing 

The  CTD  offers  both  serial  and  parallel  data  ports  for  use  in  interfacing.  In  the  CTD  data  flow 
each  byte  (8  bits)  is  preceded  by  one  bit  containing  a  logic  "zero"  and  terminated  by  two  bits,  each  of 
which  contains  a  logic  "one.”  When  a  digitized  channel  requires  two  bytes  to  convey  magnitude,  the 
first  byte  is  the  least  significant  byte.  The  original  interface,  developed  by  the  Oceanographic  Computer 
Applications  Group  in  1977,  performed  a  sample-and-hold  operation  on  an  entire  scan  of  data  holding 
the  three  16-bit  words  containing  pressure,  temperature,  and  conductivity.  These  three  words  are  then 
sent  to  the  computer  on  demand  using  three  computer  microcircuit  interface  boards,  one  for  each 
parameter.  A  driver  developed  by  the  Oceanographic  Computer  Applications  Group  for  doing  sampling 
of  successive  I/O  channels  on  HPlOOO-scries  computers  was  employed  to  control  the  I/O 
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Figure  1  —  Computer  system  aboard  USNS  Hayes  in  1977 


The  original  approach  had  some  disadvantages.  While  syncronization  problems  have  been  found 
to  be  rare,  the  "sync"  and  utility  frames  were  unavailable.  Also,  the  requirement  for  three  successive 
channels  was  inconvenient  and  occasionally  necessitated  the  use  of  an  I/O  extender  to  provide  addi¬ 
tional  computer-interface  channels.  The  extender  requires  additional  rack  space  and  creates  the  poten¬ 
tial  for  hardware  failure  aboard  ship.  Accordingly,  the  CTD  system  required  revision  to  work  with  a 
CTD  microprofiler  modified  to  include  fluorometer  measurements.  The  data  were  still  taken  at  the 
deck  unit  parallel  output  port;  however,  the  interrupt  strobes  were  "or"ed  together  and  data  was  sent  to 
the  computer  in  word  serial  mode  with  no  indication  of  which  word  (or  byte)  of  data  was  being  sent 
Thus,  the  software  had  to  determine  where  the  sync  frame  was  in  order  to  start  data  acquisition  and  to 
assure  that  the  data  flow  remained  in  sync.  The  data  were  sent  to  the  computer  through  an  I/O  board 
using  the  same  driver  as  before.  This  method  brings  in  all  the  collected  data  bytes,  but  greater  pro¬ 
gramming  care  was  required  since  the  internal  timing  of  the  CTD  within  a  scan  had  to  be  considered  in 
the  programming  to  assure  that  sufficient  time  was  available  between  words  to  perform  all  CTD  opera¬ 
tions  without  data  loss. 

The  fluorometer  CTD  system  has  the  following  data  format: 

Byte  1 :  Frame  sync.  Alternate  between  1 1 1 10000  and  00001 1 1 1 

Byte  2:  Pressure  (LSBs) 
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Byte  3:  Pressure  (MSBs) 

Byte  4:  Temperature  (LSBs) 

Byte  5:  Temperature  (MSBs) 

Byte  b:  Conductivity  (LSBs) 

Byte  7:  Conductivity  (MSBs) 

Byte  8:  Utility.  Contains  pressure  and  temperature  signs. 

Byte  9:  Fluorometer  "Signal."  All  "ones"  correspond  to  lull  scale  input. 

Byte  10:  Fluorometer  "Range."  All  "ones"  correspond  to  full  scale  input. 

As  currently  interfaced,  pressure,  temperature,  and  conductivity  are  packed  into  words  by  the 
interface  unit,  with  the  remaining  parameters  placed  in  the  lower  half  of  full  words. 

This  report  deals  with  the  CTD  software  for  the  fluorometer  CTI).  However,  both  program  sets 
are  quite  similar  and  differences  are  likely  to  be  transparent  to  a  user. 

REAL-TIME  COLLECTION  AND  DISPLAY  SOFTWARE 


The  Real-Time  Program  Set 

The  real-time  data  collection  and  display  programs  for  the  fluorometer  CTD  unit  are  described  in 
this  chapter. 


There  are  nine  real-time  programs  used.  These  are: 


CTDRT 

CTD1N 

CTDMT 

CTDP1 

CTDP2 

CTDCN 

CTDTR 

CTDMS 

CTDFL 


Obtains  data  and  schedules  other  programs. 

Initializes  parameters  prior  to  data  collection  and 
writes  an  initialization  record  on  magnetic  tape. 

Writes  data  records  on  magnetic  tape. 

Produces  real-time  plots  of  salinity  and  temperature 
vs  pressure. 

Produces  real-time  plots  of  salinity,  temperature,  and 
pressure  vs  time. 

Changes  certain  control  functions  while  data  are  being  taken. 
Writes  termination  file. 

Writes  message  to  terminal  for  CTDRT. 

Sets  fluorometer  On/Off  switch. 


System  Requirements 


The  programs  are  designed  to  operate  under  RTE-II,  RTE-III,  or  RTE-IV  (Real-Time  Executive 
operating  systems)  on  HPlOOO-series  computers  with  64K  words  of  memory.  A  background  system 
common  area  of  2000g  is  required.  The  collection  program  CTDRT  is  loaded  as  a  foreground  program 
using  reverse  (background)  common.  The  remaining  programs  are  loaded  as  background  programs 
using  background  common. 


The  programs  listed  in  this  section  utilize  the  HP2647  terminal  and  HP2631G  printer.  The 
HP2631G  is  also  used  to  produce  hard  copy  of  terminal  plots.  The  programs  require  minor  revisions  if 
other  peripherals  are  used  since  the  means  for  producing  hard  copy  of  plots  is  device  dependent  and 
also  the  screen  size  of  the  HP2647  requires  a  minor  scaling  modification  of  the  plot. 
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Data  File  Structure 

One  CTD  lowering  (or  towing)  results  in  a  data  tape  file  with  the  following  format. 

I  INITIALIZATION  RECORD  |  EOF  I  DATA  RECORDS  I  EOF  |  TERM  RECORD  I  EOF  I  EOF  | 

Prior  to  the  next  event,  the  program  F2EOF  is  run  to  position  the  tape  between  the  double  end-of-file 
marks  and  the  next  initialization  record  is  written  over  the  last  EOF.  Hence  a  tape  with  several  CTD 
events  will  have  a  double  EOF  only  after  the  last  event.  A  sequence  of  three  events  would  appear  as 


E 

E 

E 

E 

E 

E 

E 

F. 

E 

E 

lnit. 

0 

Data 

0 

Term. 

0 

lnit. 

0 

Data 

0 

Term. 

0 

lnit. 

0 

Data 

0 

Term. 

0 

0 

Record 

F 

F 

Record 

F 

Record 

F 

F 

Record 

F 

Record 

F 

F 

Record 

F 

F 

The  initialization  record  is  a  47-word  header  record  containing  identification  parameters  for  the 
event: 


Words 

1-6: 

12  Character  Identifying  ASCII  String 

Word 

7: 

Station  Number 

Word 

8: 

Year 

Word 

9: 

Month 

Word 

10: 

Day 

Word 

11: 

Lowering  Number 

Word 

12: 

Instrument  Number 

Words 

13-24: 

Conversion  Constants.  Cl,  . . . ,  C6 

Words 

25-30: 

Pressure  Axis  Constants 

Words 

31-36: 

Salinity  Axis  Constants 

Words 

37-42: 

Temperature  Axis  Constants 

Words 

43-44: 

Latitude 

Words 

45-46: 

Longitude 

Word 

47: 

Zero 

The  initialization  record  is  followed  by  an  EOF  and  then  the  CTD  data.  Each  data  record  contains 
456  words  (approximately  2  seconds  of  data)  consisting  of  8  words  of  header  information  followed  by 
448  data  words: 


Word 

Word 

Word 

Word 

Word 

Word 

Word 

Word 

Words 


1:  Record  Number 

2:  Time  (Hours  =  24  x  Days  +  Hours) 

3:  Time  (Seconds  =  60  x  Minute  +  Seconds) 

4:  Time  (Tens  of  milliseconds) 

5:  Special  Segment  ID  Marker 

6:  Fluorometer  on/off  switch 

7:  0 

8:  0 

9-456:  Digital  Data 


The  last  data  record  is  followed  by  an  EOF  and  an  8-word  termination  record: 


Word 

1: 

-1 

Word 

2: 

Time  Hours 

Word 

3: 

Time  Seconds 

Word 

4: 

Termination  Method  Indicator 

Words 

5-8: 

0 

Finally,  a  double  EOF  is  written  on  the  tape. 
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Data  Acquisition  and  Display 

The  da.  i-acquisition  programs  .ire  started  by  running  the  program  CTDRT  (a  full  description  of 
questions  and  tnswers  is  given  in  the  next  chapter)  CTDRT  then  schedules  programs  CTDMS  for 
execution.  CTDMS  asks  the  operator  whether  the  operation  will  be  a  CTD  lowering  or  tow  and  if  tem¬ 
perature  or  potential  temperature  is  to  be  a  plotting  parameter  For  a  lowering,  temperature  (or  poten¬ 
tial  temperature)  and  salinity  are  both  plotted  against  pressure,  while  in  the  case  of  a  tow,  temperature 
(or  potential  temperature),  salinity,  and  pressure  are  plotted  as  a  function  of  time.  Note  that  pressure 
in  decibars  is  very  nearly  identical  to  depth  in  meters. 

Now  that  the  plot  type  has  been  defined,  program  CTDRT  schedules  programs  CTDIN  for  execu¬ 
tion  to  obtain  identifying  parameters  and  plot  scale  parameters.  These  parameters  are  written  to  mag¬ 
netic  tape  in  a  47-word  header  record  for  the  CTD  cast.  The  calibration  constants  are  used  to  convert 
the  CTD  measurement  into  pressure,  temperature,  and  salinity.  If  Mt.  M2.  A/y  are  the  real  values  of 
the  measurements  obtained  from  the  CTD  unit  for  conductivity  pressure  and  temperature,  then 

Pressure  =  Cl  *  (M2  —  C 2) 

Temperature  =  C 3  +  (C4  *  M}) 

Conductivity  =  C5  *  A/,  *  (1  +  C 6  *  Temperature) 

Cl  is  the  only  constant  with  which  the  user  need  be  concerned;  it  depends  upon  the  pressure  sen¬ 
sor  used  within  the  CTD.  The  axes  parameters  determine  the  parameter  limits  for  the  real-time  CTD 
display.  Points  outside  the  user-defined  axes  limits  will  not  appear  on  the  real-time  plot,  so  the  user 
should  exercise  care  in  assuring  that  these  limits  span  the  data  without  being  so  wide  that  the  display 
loses  meaning.  Note,  however,  that  these  limits  affect  only  the  real-time  display  and  not  the  collection 
of  the  data  or  the  process  of  writing  the  data  records  onto  magnetic  tape.  Thus,  if  the  date  is  outside 
the  user-defined  limits,  the  ability  to  view  the  data  in  real-time  is  lost  but  the  collection  process  will  not 
be  affected. 

Once  the  initialization  parameters  have  been  obtained,  the  program  CTDIN  displays  the  values 
and  asks  if  the  values  are  correct.  If  a  mistake  has  been  made  while  inputting  the  parameters,  the  user 
answers  NO  and  the  questions  are  repeated.  If  they  are  correct,  CTDIN  writes  the  initialization  record 
and  an  end-of-file  (EOF)  onto  magnetic  tape  and  terminates,  which  allows  the  scheduling  program 
CTDRT  to  resume. 

Depending  upon  whether  a  lowering  or  towing  real-time  display  has  been  requested,  CTDRT  next 
schedules  either  program  CTDP1  or  CTDP2.  In  either  case  the  appropriate  axes  are  drawn  on  a  plot¬ 
ting  terminal,  with  the  axes  parameters  obtained  by  CTDIN.  When  these  axes  have  been  drawn, 
CTDP1  (or  CTDP2)  suspends  execution  with  the  point  of  suspension  being  retained.  This  enables 
plotting  to  occur  without  the  screen  being  erased  each  time  a  pair  of  points  are  to  be  drawn.  Program 
CTDRT  now  has  program  CTDMS  print  the  message  "PAUSE— TYPE:  ‘GO,  CTDRT'  FOR  DATA 
ACQUISITION"  on  the  system  console  and  then  suspends.  This  completes  the  initialization  phase. 
The  operator  then  types  "GO,  CTDRT'  on  the  system  console  to  initiate  data  collection. 

Figures  2  and  3  were  made  aboard  the  USNS  Lynch  in  October  1977  near  Bermuda  and  show  typi¬ 
cal  plots  for  the  lowering  mode.  For  Fig.  2  the  CTD  was  lowered  near  the  surface  while  in  Fig.  3  the 
device  "yo-yoed”  near  the  bottom.  Figures  4  and  5  were  made  during  the  return  run  from  Bermuda  to 
Newport,  Rhode  Island,  and  show  the  real-time  display  for  the  towing  mode  Here  pressure,  tempera¬ 
ture,  and  salinity  are  displayed  as  a  function  of  time. 
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The  basic  programming  technique  used  for  data  acquisition  and  processing  is  double  buffering. 
The  data  are  acquired  in  (approximately)  2-second  buffers  and  the  processing  of  the  previous  buffer 
must  be  accomplished  within  the  2  seconds  allotted  for  obtaining  the  new  data.  Since  the  CTD  scan 
rate  is  31.25  hertz,  one  scan  takes  32  milliseconds.  As  previously  discussed,  the  CTD  interface  sends 
the  data  words  in  a  serial  format.  The  CTD  interface,  in  one  scan,  sends  the  computer: 

Word  1:  Sync  Frame 

Word  2:  Pressure 

Word  3:  Temperature 

Word  4:  Conductivity 

Word  5:  Utility 

Word  6:  Fluorometer  ''Signal’' 

Word  7:  Fluorometer  "Range” 

where  those  words  which  contain  only  one  byte  of  information  have  the  data  in  the  LSBs. 

Returning  to  the  data  processing  flow,  program  CTDRT  compares  words  against  the  CTD  sync- 
frame  values  until  the  sync  frame  is  recognized  It  next  determines  which  of  the  two  buffers  is  to  be 
filled.  Once  this  is  done,  program  CTDRT  schedules  program  CTDMT  and  requests  448  words  of  data 
from  the  CTD  interface.  Thus,  while  the  foreground  program  CTDRT  obtains  the  next  2  seconds  of 
data,  the  background  programs  (CTDMT  and  either  CTDP1  or  CTDP2)  processes  the  previous  2 
seconds  of  data.  CTDMT  keeps  track  of  the  buffer  times  and  fills  in  the  appropriate  entries  in  the 
header— the  first  8  words  of  the  buffer  as  given  above.  The  record  is  then  dumped  onto  magnetic  tape, 
and  every  fourth  buffer  the  plotting  program  CTDP1  (or  CTDP2)  is  scheduled  to  plot  a  point  and  print 
out  on  the  line  printer  and  system  console  the  values  of  temperature,  pressure,  and  salinity.  These 
values  are  printed  even  if  a  point  is  not  plotted  due  to  one  of  the  values  being  off  the  plot  scale.  When 
program  CTDP1  (or  CTDP2)  is  scheduled  for  execution,  it  takes  the  first  three  data  points  in  the 
buffer  and  uses  the  constants  Cl.  C2,  ....  C 6  to  obtain  pressure,  temperature,  and  conductivity. 
Salinity  is  computed  using  an  algorithm  of  Fofonoff.  These  points  are  then  scaled  for  plotting  and,  if 
within  the  user  defined  plot  axes  limits,  they  are  plotted  on  the  real-time  plot  display.  The  values  are 
printed  out  as  mentioned  above.  If  the  CTD  is  being  operated  in  the  tow  mode  each  plot  will  contain 
1/2  hour  of  data.  At  the  end  of  the  half-hour  program  CTDP2  will  automatically  copy  the  plot  onto  a 
hard-copy  unit  (the  263 1G  terminal,  a  Tektronix  hardcopier,  etc.,  depending  upon  the  system  in  use). 
The  axes  are  then  redrawn  and  reannotated  with  the  same  value  as  previously  used  and  the  plot  contin¬ 
ues. 

The  remaining  programs  are  used  for  various  control  functions.  Program  CTDFL  sets  the 
fluorometer  on/off  switch  (Word  6  of  the  buffer),  it  should  be  noted  that  this  switch  is  for  conveni¬ 
ence  in  locating  fluorometer  data,  but  is  not  used  to  alter  the  data,  hence,  operator  failure  to  properly 
set  the  switch  will  not  affect  the  data  tape.  The  program  CTDCN  is  used  for  several  purposes.  It  can 
be  used  to  put  a  marker  in  word  5  of  the  data  record  in  order  to  identify  special  segments  of  data  for 
easy  future  location.  Alternately,  program  CTD  is  used  to  put  the  CTD  programs  in  "idle"  mode  and  to 
terminate  data  collection.  Occasionally,  in  the  course  of  a  CTD  lowering,  events  such  as  cable  snags 
will  necessitate  holding  the  CTD  at  a  fixed  depth  for  a  certain  period,  and  it  has  been  found  useful  to 
have  the  capability  of  suspending  data  collection  while  such  problems  are  being  corrected.  Accordingly, 
one  mode  of  operation  of  program  CTDCN  will  put  all  programs  into  a  suspend  state  in  such  a  way  that 
the  operator  can  resume  data  collection  by  typing  a  command  on  the  system  console.  CTDCN  is  also 
used  to  terminate  data  collection.  Full  details  of  how  to  use  program  CTDCN  are  given  the  next  sec¬ 
tion 

When  program  CTDCN  is  executed  in  either  its  termination  or  its  abort  mode,  it  assures  that  the 
real-time  programs  are  properly  terminated,  writes  an  F.OF  onto  the  magnetic  tape,  and  then  schedules 
program  C  TDTR  which  writes  an  8-word  termination  file  and  two  HOI  s.  This  completes  the  data  col¬ 
lection  process  for  one  event  For  further  collection  using  the  same  magnetic  tape,  the  tape  is  posi¬ 
tioned  between  the  EOF's  and  the  process  repeated. 


8 


NRL  REPORT  8444 


Real-Time  Programs— FORTRAN  Source  Listings 


0001  FTN4.L 

0002  PROGRAM  CTDRT( 10, 15) 

0003  C 

0004  C  REAL  TIME  PROGRAMS  FOR  COLLECTION  OF  CTD  DATA  WITH  FLOROMETER 

0005  C  CAPABILITY.  DOUBLE  BUFFERED  DATA  TAKEN  EVERY  TUO  SECONDS. 

0006  C 

0007  C  REAL  TIME  PROGRAMS  USED  ARE: 

0008  C  CTDRT  -  SCHEDULES  OTHER  PROGRAMS  AND  CONTROL  PARAMETERS 

0009  C  CTDIN-  INITIALIZES  PARAMETERS  AND  WRITES  IN  IT  FILE  ON  MT 

0010  C  CTDMT  -  WRITES  DATA  BUFFERS  TO  MT 

0011  C  CTDP1-  REAL  TIME  PLOT  OF  SAL,  TEMP  VS.  PRESSURE 

0012  C  CTDP2-  REAL  TIME  PLOT  OF  SAL.  TEMP.  PRESS  VS.  TIME 

0013  C  CTDTR-  WRITES  TERMINATION  FILE 

0014  C  CTDCN-  USES  COMMON  PARAM  ISTOP  TO  CONTROL  FUNCTIONS 

0015  C  CTDMS-  USED  TO  WRITE  MESSAGES  FOR  CTDRT 

0016  C  CTDFL-  SETS  FLOROMETER  ON/OFF  SWITCH 

0017  C  OFF  LINE  PROGRAMS  USED  ARE: 

0018  C  PLTBK  -  USED  FOR  POST  PROCESSING  OF  DATA 

0019  C  F2E0F  -  POSITIONS  MT  BETWEEN  DBL  EOF  FOR  NEW  RT  DATA 

0020  C  MTPOS  -  FINDS  STA  &  LOU  NUMBERS  AND  POSITIONS  MT  AT  DATA 

0021  C 

0022  C  DATA  FILE  FORMAT 

0023  C  WORD  1  RECORD  NUMBER 

0024  C  WORD  2  TIME (HOURS) 

0025  C  WORD  3  TIME (SECONDS) 

0026  C  WORD  4  TIME ( 10XMSEC) 

0027  C  WORD  5  ID  MARKER  FOR  SPECIAL  SEGMENTS  -  SEE  CNTL 

0020  C  WORD  6  FLOROMETER  ON/OFF  0-OFF  1-ON 

0029  C  WORDS  7-8  0 

0030  C  WORDS  9-456  DIGITAL  DATA 

0031  C 

0032  C  SEE  CTDIN  FOR  47  WORD  INITIALIZATION  FILE  AND  CTDTR  FOR  8  WORD 

0033  C  TERMINATION  FILE 

0034  C 

0035  C  SCIENTIFIC  SUBROUTINES  USED  ARE  TAKEN  FROM  WOODS  HOLE  CTD 

0036  C  PROGRAM  AND  ARE  LISTED  WITH  THE  USER  PROGRAM 

0037  C 

0038  C  PROGRAMMER:  LARRY  ROSENBLUM 

0039  C  VERSION:  SEPT  20,1979 

0040  C 

0041  DIMENSION  NAME (3) , NAME  1 (3) .NAME2 (3) ,NAME3(3) , ITIME (5) 

0042  DIMENSION  NAME4(3) , ISYNC (7) .NAME5 (3) , IPRAM(5) 

0043  COMMON  IREC, ISTOP, ISCAN.NSCAN, IFLOR 

0044  COMMON  IBUFK456),  IBUF2C456) 

0045  COMMON  C1,C2,C3,C4,C5,C6, INSTR 

0046  COMMON  ISHIP(6) , 1ST, IYR, I MO, I DAY, LOU 

0047  COMMON  PAXL.P IV,  P INCR, TAXL, TIV, TINCR, SAXL, S I V, S INCR 

0048  DATA  NAME72HCT.2HDM.2HT  7 

0049  DATA  NAME172HCT.2HDI.2HN  7 

0050  DATA  NAME272HCT.2HDT.2HR  7 

0051  DATA  NAME372HCT.2HDP.2Ht  7 

0052  DATA  NAME4/2HCT, 2HDP, 2H2  7 
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3053 

0054 

0055 

0056 

0057 

0050 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 


DATA  NAhE5x2HCT.2HDM.2HS  / 

LUCTD-10 
MSK 1-1774008 
MSK2-  3778 

C  CHOOSE  TYPE  OF  PLOT  WANTED 

CALL  EXEC (9, NAMES. 1) 

CALL  RMPARC IPRAM) 

IANS-IPRAMC2) 

IANS1-IPRAMC3) 

I  STOP— 1 

IREC-1 

IFLOR-0 

C  SCHEDULE  IN  IT  AND  PLT1  OR  PLT2 

CALL  EXEC (9, NAME  1. IANS) 

IF(  IANS.EQ. 1 ) CALL  EXECC 10. NAME3, IANS1) 

IF ( IANS . F  2) CALL  EXECC 10, NAME4,  IANS1) 

C  WAIT  U  iL  AXES  HAVE  BEEN  DRAUN  &  EXPT  READY  TO  COLLECT  DATA 

CALL  EXEC (9, NAMES. 2) 

PAUSE 

C  GET  INTO  SYNC  WITH  THE  SYNC  FRAME 

2  CALL  EXECC l.LUCTD, ISYNC. 1) 

I TST- I AND (I SYNC <1 ) . MSK 1 > 

IFC ITST.NE. 170000B)GO  TO  2 
CALL  EXECC l.LUCTD, ISYNC,  6) 

C  DETERMINE  WHICH  BUFFER  IS  TO  BE  FILLED 

12  KPAR-IREC-2*(IRECX2) 

IF (KPAR.EO.0)GO  TO  48 
C  ODD  RECORD  NUMBER 

IBUF1CD-IREC 
IBUF1(6)-IFL0R 
C  SCHEDULE  MAGTP 

IF C  IREC . NE .  1 ) CALL  EXEC C 18, NAME. IANS, IANS 1 ) 
r  QOIIIRF  DATA 

60  CALL  EXECC l.LUCTD, IBUF1 (9) ,448) 

IREC* IREC+1 

IF  ( ISTOP.GE.0)  GO  TO  99 
IFCIREC.LT. 4)G0  TO  12 
ITST1-IANDCMSK1. IBUF2C9) ) 

IFC ITST1 .EQ. 170000B) GO  TO  12 
IFC ITST1 .EQ.7400B)GO  TO  12 
ITST1-IANDCMSK2.  IBUF2C9)) 

IFC 1 • ~Tl . EO.360B)GO  TO  12 
IFCITST1.EO. 17B)G0  TO  12 
GO  TO  2 

C  EVEN  RECORD  NUMBER 

40  IBUF2  Cl)*  IREC 
IBUF2C6) *IFLOR 
C  SCHEDULE  MAGTP 

CALL  EXECC 10, NAME, IANS, IANS1) 

C  AQUIRE  DATA 

CALL  EXECC l.LUCTD. IBUF2C 9). 448) 

IREC- IREC+1 

IFC ISTOP.GE.0)  GO  TO  99 
IFCIREC.LT. 4)GO  TO  12 
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010?  ITST1-IAND(MSK1.  IBUFK9)) 

0108  1FUTSTI.EQ.  1?0000B)GO  TO  12 

0109  I F  C ITST1 . EQ . 00?400B) GO  TO  12 

0110  ITST1-IAND(MSK2, IBUF 1(9)) 

0111  IF < ITST1 .EQ.360B)GO  TO  12 

0112  IF( ITST1 . EQ. l?B)GO  TO  12 

0113  GO  TO  2 

0114  C  END  OR  IDLE 

0115  99  IF C  ISTOP. NE. 2)  GO  TO  90 

0116  CALL  EXEC(9,NAME5.3) 

0117  PAUSE 

0116  ISTOP  — 1 

0119  GO  TO  2 

0120  98  CALL  EXEC(9.NAME2, IANS) 

0121  END 

0122  ENDS 


0001  FTN4.L 

0002  PROGRAM  CTDIN 

0003  C 

0004  C  PROG  IS  CALLED  (WITH  WAIT)  BY  CTDRT  TO  INITIALIZE  PARAMS. 

0005  C  PARAMETER  NAMES  AND  INITIALIZATION  FILE  WERE  MADE  TO 

0006  C  CORESPOND  TO  THE  UOODS  HOLE  PROGRAM.  SOME  PARAMETERS 

000?  C  ARE  NOT  USED  IN  PROGRAM.  4?  WORD  FILE  IS  WRITTEN  TO  MT 

000B  C  ON  COMPLETION. 

0009  C  INSTR- INSTRUMENT  NUMBER 

0010  C  C1-C6  ARE  CONSTS  USED  TO  CONVERT  MEAS.  TO  PRES. TEMP,  a  COND 

0011  C  AXES  PARAMS  ARE  AX  LENGTH, IN  IT.  VALUE,  a  INCR/INCH.  AXES 

0012  C  LENGHTH  IS  6"  FOR  X  AXIS  AND  5"  FOR  Y  AXIS. 

0013  C 

0014  C  FILE  FORMAT: 

0015  C  WORDS  1-6  12  CHARACTER  INPUT  I.  D.  STRING 

0016  C  WORD  ?  STATION  NUMBER 

0017  C  WORDS  B- 10  DATE:  YEAR,  MONTH.  DAY 

0018  C  WORDS  11-12  LOWERING  *,  INSTRUMENT  # 

0019  C  WORDS  13-24  CALIBRATION  CONSTS  (REAL  VALUED) 

0020  C  WORDS  25-30  PRESSURE  AXIS  CONSTS  (REAL  VALUED) 

0021  C  WORDS  31-36  TEMPERATURE  AXIS  CONSTS  (REAL  VALUED) 

0022  C  WORDS  37-42  SALINITY  AXIS  CONSTS  (REAL  VALUED) 

0023  C  WORDS  43-46  LATITUDE  a  LONGITUDE  (REAL  VALUED) 

0024  C  WORD  4?  0 

0025  C 

0026  C  PROGRAt't'ER :  LARRY  ROSENBLUM 

0027  C 

0028  DIMENSION  IBUF5(47) , IPRAM(5) 

0029  COMMON  IREC. ISTOP, ISCAN.NSCAN, 1FLOR 

0030  COMMON  IBUF 1 (456) . IBUF2 (456) 

0031  COMMON  Cl. C2,C3.C4,C5,C6, INSTR 

0032  COMMON  ISHIP(6). 1ST. IYR. IMO, IDAY.LOW 

0033  COMMON  PAXL.PIV.PINCR.TAXL.TIV.TINCR.SAXL.SIV.SINCR 

0034  EQUIVALENCE  (XC 1 . IBUF5( 13) ) , (XC2, IBUF5( 15) ) 
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0035 

EQUIVALENCE  CXC3. IBUF5C 17) ) , (XC4, IBUF5 ( 19) ) , CXC5, IBUF5C21)) 

0036 

EQUIVALENCE  (XC6. IBUF5<23> ) , (XPAXL, IBUF5C25) ) , CXPIV, IBUF5C27)) 

0037 

EQUIVALENCE  (XPINCR, IBUF5 C29) ) , (XTAXL. I8UF5C31) ) , (XTIV, IBUF5C33)) 

0038 

EQUIVALENCE  (XTINCR, IBUF5 (35) ) , CXSAXL, IBUF5C3?)),  (XSIV,  IBUF5C39)) 

0039 

EQUIVALENCE  (XSINCR,  IBUF5C41) ) ,  (XL  AT,  IBUF5C43) ) , (XLON,  IBUF5C45)) 

0040 

CALL  RMPARC IPRAM) 

0041 

C 

DEFAULT  VALUES 

0042 

5 

INSTR-1 

0043 

Cl-0. 1 

0044 

C2-0.0 

0045 

C3-0.0 

0046 

C4-.0005 

0047 

C5-0.001B 

0048 

C6-0.0 

0049 

PAXL-5. 

0050 

PIV-0. 

0051 

PINCR-1000. 

0052 

TAXL-6. 

0053 

IF ( IPRAMC 1) . EQ . 2) TAXL-5 . 

0054 

TIV-1. 

0055 

TINCR-3.0 

0056 

IF ( IPRAMC  1) .EQ.2)TINCR-4.0 

0057 

SAXL-6. 

0058 

IF ( IPRAMC 1) . EQ.2)SAXL-5. 

0059 

SIV-32. 

0060 

IF( IPRAMC 1) . EQ.2)S IV-33. 

0061 

SINCR-1.0 

0062 

IBLNK-20040B 

0063 

DO  37  1-1.6 

0064 

ISHIPC I) -I8LNK 

0065 

37 

CONTINUE 

0066 

URITECl. 10) 

0067 

10 

FORMAT ( "PROGRAM  CTDIN  RUNNING") 

0068 

C 

READ  IN  PARAMETER  VALUES 

0069 

URITEC 1, 12) 

0070 

12 

FORMAT ( "ENTERS  SHIP  +  CRUISE  *  (12  CHARACTERS)") 

0071 

READU.  13)  (ISHIPC  I).  1-1,6) 

0072 

13 

F0RMAT(6A2) 

0073 

URITE(1. 14) 

0074 

14 

FORMAT! "ENTER  DATE (YR. MO. DAY) ") 

0075 

READ ( 1 , *) IYR, IMO,  IDAY 

0076 

URITECl. 16) 

0077 

16 

FORMAT! "ENTER  STATION  NUMBER  AND  LOUERING  NUMBER") 

0078 

READ ( 1 ,  #) 1ST.  LOU 

0079 

URITECl. 17) 

0080 

17 

FORMAT ( "ENTER  LATITUDE:  DEG,  MIN  &  FRACT.  OF  MIN  XXXX. , XX. XXX") 

0081 

READ(1.*)XLAT1,XLAT2 

0082 

URITECl. 19) 

0083 

19 

FORMATC "ENTER  LONGITUDE: DEG.  MIN  8.  FRACT.  OF  MIN  XXXX. .XX. XXX") 

0084 

READ ( 1 , *) XLON 1 . XL0N2 

0085 

XLAT-60.0*XLAT1+XLAT2 

0086 

XLON -60 . 0*XLON 1 +H.0N2 

0087 

URITECl, IB) 

0088 

18 

FORMATC "USE  DEFAULT  VALUES,  1-YES  2-NO") 

0089 

READ! l,*)KANS 

0090 

IFCKANS.EQ. l)GO  TO  99 

0091 

URITECl. 20) 

0092 

20 

FORMATC "DEFAULT  VALUES  GIVEN  IN  PARENTHS") 

0093 

URITECl, 22) 
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0094 

0095 

0096 

009? 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

010? 

0108 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

014? 

0148 

0149 

0150 

0151 

0152 


22  FORMAT! ■ INSTRUMENT  •  (1) ■) 

READ! 1,*) INSTR 
UR1TE! 1,26) 

26  FORMAT! "CAL  CONST  PRESSURE  C1.C2  !0. 1,0.0)*) 

READ! 1,*)C1.C2 
UR  I TEC  1,28) 

28  FORMAT ! "CAL  CONST  TEMP  C3.C4  10.0,0.0005) *) 

READ! 1,*)C3,C4 
WRITE! 1,30) 

30  FORMAT! "CAL  CONST  COND  C5,C6  (0.001,0.0)’) 

READ! 1,*)C5,C6 
WRITE! 1,31) 

31  FORMAT! ’FOR  EACH  AXIS  THE  AXIS  LENGTH  IS  GIVEN  IN  PARENTHESIS’) 
WRITE! 1,33) 

33  FORMAT! "ENTER  THE  INITIAL  VALUE  AND  THE  INCREMENT/INCH’) 

WRITE! 1, 32) PAXL 

32  FORMAT! "PRESSURE(AX  LNTH-’, 12, ’)ENTERs  IN  IT  VALUE.  INCR  INCH’) 
READ! 1,*)PIV, PINCR 
WRITE! 1. 34) TAXL 

34  FORMAT ( ’TEMP! AX  LNTH-’. I 2,’) ENTER  INIT  VALUE.  INCR/INCH’) 

READ(1,*)TIV.TINCR 
WRITE! 1. 36) SAXL 

36  FORMAT (’SAL (AX  LNTH-’, I 2,’) ENTER  INIT  VALUE.  INCR. INCH*) 

READ! 1,*)SIV,SINCR 
C  WRITE  OUT  VALUES 

99  DO  98  1-1,6, 5 

WRITE! 1,60) (ISHIP(J).J-1,6> 

60  FORMAT! IX, 6A2) 

WRITE! 1,62)  IYS, IMO, IDAY, 1ST, LOU 
62  FORMAT!  IX,  ’DATE  ’,314,’  STA  M4,  ’  LOU  ’,14) 

WRITE! 1,64) INSTR 

64  FORMAT! IX, ’INSTR  *,  12) 

WRITE! I,66)C1,C2,C3,C4, C5.C6 

66  FORMAT! IX, ’CAL  CONSTS  *,6F9.4) 

URITE(I.68)PAXL,PIV,PINCR,TAXL,TIV,TINCR.SAXL,SIV.SINCR 
68  FORMAT! IX, ’AXES  CONSTS  ’.9F6.2) 

UR  I TE ( 1 , 70) XLAT 1 . XLAT2. XLON 1 , XL0N2 
70  FORMATCLAT:  ",F4.0,  ’DEG",F6.2,  ’MIN’.4X.  "LON:  \F4.0,  ’DEG’, 

1  F6.2. "MIN’) 

98  CONTINUE 

WRITE! t. 105) 

105  FORMAT! "ARE  INPUT  PARAMS  O.K.:  1-VES,  2-NO.  3-YES  -  NO  MT  FILE’) 
READU.*)  IANS 
IF! IANS.EQ.2)G0  TO  5 
C  WRITE  FILE  OUT  TO  TAPE 

IF ( IANS.EQ.3)  GO  TO  199 
DO  110  1-1,6 
IBUF5! I) -ISHIP!  I) 

110  CONTINUE 

IBUF5 (7) ■ 1ST 
IBUF5 (8) ■ IYR 
!BUF5(9) * IMO 
IBUF5! 10) -IDAY 
IBUF5 (11) -LOU 
IBUF5  C 12) * INSTR 
IBUF5 (47 ) -0 
XC 1 -C  1 
XC2-C2 
XC3-C3 
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0153 

0154 

0155 

0156 

015? 

0158 

0159 

0160 

0161 

0162 

0163 

0164 

0165 

0166 

016? 

0168 

0169 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

001? 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

8029 

0030 

0031 


XC4-C4 

XC5-C5 

XC6-C6 

XPAXL-PAXL 

XP1V-PIV 

XPINCR-PINCR 

XTAXL-TAXL 

XTIV-TIV 

XTINCR-TINCR 

XSAXL-SAXL 

XSIV-SIV 

XSINCR-SINCR 

C  URITE  RECORD  AND  EOF 

CALL  EXEC (2, 8, IBUF5, 47) 
CALL  EXEC(3,8+100B) 

199  END 
END* 


FTN4.L 

PROGRAM  CTDMT(3,90) 

C 

C  PROGRAM  TO  URITE  OUT  DOUBLE  BUFFERED  DATA  COLLECTED  BY 

C  CTDRT  TO  MAG.  TAPE.  ALSO  SCHEDULES  RT  PLOTTING  PROG  AND 

C  GETS  TIME. 

C 

C  PROGRAMMER:  LARRY  ROSENBLUM 

C  VERSION  SEPT  22  1979 

C 

C  NOTE:  1ST  RECORD  MAY  HAVE  BAD  TIME 

C 

DIMENSION  ITIMEC5),  IPRAMC5) .  NAME3C3) .NAME4C3) 

COMMON  IREC, ISTOP, ISCAN, NSCAN, IFLOR 
COMMON  IBUF 1(456). IBUF2C456) 

COMMON  C1.C2.C3.C4.C5.C6. INSTR 
COMMON  ISHIP<6>. 1ST,  IYR,  IMO, IDAY.LOU 
COMMON  PAXL.PIV, PINCR.TAXL.TIV.TINCR.SAXL.SIV, SINCR 
DATA  NAME3/2HCT, 2HDP, 2H 1  / 

DATA  NAME4/2HCT.2HDP.2H2  / 

5  CALL  RMPARC IPRAM) 

IANS  -IPRAM(l) 

IANS1-IPRAM(2) 

KPAR» IREC-2*( IREC/2) 

C  GET  CURRENT  TIME 

CALL  EXECdl.  ITIME) 

IHR*24#ITIME(5)+ITIME(4) 

ISEC*60*ITIME (3) +ITIME (2) 

IMSC-ITIME(l) 

C  SCHED.  PLOTTING  PROG  EVERY  4TH  TIME 

TEST*  IREC-4>K(  IREC/4) 
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0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0006 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0016 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 


IF ( TEST. NE. 3) GO  TO  15 
IF C IANS.EQ. 1)CALL  EXECC 10.NAME3. IANS1) 
IF(  IANS.EQ. 2) CALL  EXECU0.NAME4, IANS1) 
15  IF  (KPAR.EQ.0)  GO  TO  10 
I8UF2C2) -IHR1 
IBUF2C3) *ISEC1 
IBUF2(4) • IMSC 1 
CALL  EXEC(2,8, IBUF2.456) 

GO  TO  20 

10  IBUF 1 (2) •  IHR1 

IBUF1 (3) "ISEC1 
IBUF 1(4)- IMSC 1 
CALL  EXEC(2,8, IBUF1.456) 

C  SAVE  OLD  TIME 

20  IHR1-IHR 

ISEC1-ISEC 
IMSC 1- IMSC 

IF ( ISTOP.EO. 1 .OR. ISTOP. EQ.3)G0  TO  99 
CALL  EXEC (6. 0.1) 

GO  TO  5 
99  END 

END* 


FTN4.L 

PROGRAM  CTDP1 
C 

C  REAL  TlfE  PLOTTING  PROGRAM  FOR  CTD. 

C  USES  VALUE  OF  IREC  FROM  COMMON  TO  PLOT  ONE  POINT  EVERY  8TH 

C  BUFFER.  USES  VALUE  OF  ISTDP  TO  END  OR  REINITIALIZE. 

C  RELOCATABLES:  XLRSAL,  XR,  XPOTMP 

C 

C  PROGRAMMER:  LARRY  ROSENBLUM 

C  VERSION:  SEPT  18.  1979 

C 

D I MENS  ION  DAT ( 3) . NAMES ( 4) . NAMET  C 6 ) , NAMEP ( 4) . NASTA ( 2 ) , NALOUC 2 ) 
DIMENSION  IPRAMC5) ,NAP0T(7) 

COMMON  IREC. ISTOP. ISCAN.NSCAN. IFLOR 
COMMON  IBUF 1 (456), IBUF2 (456) 

COMMON  C1,C2.C3.C4,C5,C6, 1NSTR 

COMMON  ISHIP(6) . 1ST, IYR, I MO, I DAY. LOU 

COMMON  PAXL . P I V. P INCR , TAXL , TIV, TINCR, SAXL . S I V, S I NCR 

DATA  NAMEP-'2HPR,2HES, 2WSU.2HRE/ 

DATA  NAMES/2HSA.2HLI.2HNI.2HTY/ 

DATA  NAMET^2HTE, 2HMP, 2HER, 2HAT.2HUR. 2HE  ✓ 

DATA  NAP0T/2HP0.2HTE, 2HNT.2HIA.2HL  .2HTE.2HMP/ 

DATA  NASTAY2HST.2HA  / 

DATA  NALOU/2HL0.2HU  / 

CALL  RMPAR(IPRAM) 
lANSl-IPRAM(l) 

CALL  EXEC (22,0) 
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0028  7 

0029  C 

0030 

0031 

0032 

0033  C 

0034 

0035 

0036 

003? 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

004? 

0048 
0049 
0050 
0051 
0052 
0053 
0054  C 
0055 
0056  C 
005?  30 

0058  C 
0059 
0060 
0061 
0062 
0063 
0064 
0065 
0066 
006?  C 
0068  C 
0069  C 
00?0  C 
00?  1  C 
0072  C 
00?3 
0074 

0075  7?9 

0076 
007?  C 
0078  C 
0079 
0080 
0081  21 
0082 

0083  98 

0084 

0085  99 

0006 


CALL  PLTLU(?) 

DRAU  AXES.  LABEL.  AND  ANNOTATE 
NSTIC” 1NT (SAXL) +1 
NPTIC* 1 NT (PAXL) +1 
NTTIG*  1NT  (TAXL) +1 

NEXT  LINE  TO  FUDGE  FOR  2647  TERMINAL 
PAXL-0.8*PAXL 

CALL  AX(0.6,0.4, SAXL,NSTIC,3,0.0) 

CALL  AX(0 . 6, PAXL+0 . 4. PAXL, NPTIC, 3, -90 . ) 

CALL  AXC.6+TAXL, .4+PAXL,TAXL,NTTIC,3,-180.) 

CALL  LABEL (0.6, 0.2, . 1,SAXL,NSTIC,SIV.SINCR,2.0.0) 

CALL  SYM8(3.2,0.0, . 1 , NAMES, 0 . 0,  8) 

CALL  LABEL (0.3, PAXL +0.4, . 1 ,-PAXL,  NPTIC, P IV.P INCR.2,  +90. ) 

CALL  SYM8(0. 1,2.4, . l.NAMEP, 90. 0,8) 

CALL  LABEL( .6, . 6+PAXL, . 1 , TAXL.NTTIC, TIV, TINCR,  2,  0. ) 

IF ( IANS  1 . EQ .  1 ) CALL  SYM8(3.00, .B+PAXL,  . 1.NAMET.0.0,  12) 

I F ( I ANS  1 . EQ . 2 ) CALL  SYMB (3.00,. 8+PAXL, . 1 . NAPOT, 0.0,14) 

CALL  SYMB ( . 7+S AXL , - . 4+P AXL , . 1 , NASTA , 0 . 0 , 4) 

CALL  SYMB ( . 7+SAXL, - . 7+PAXL,  . l.NALOU,0.0.4) 

RIST-R( 1ST) 

RLOLJ-R(LOU) 

CALL  LABEL(1.2+SAXL.-.4+PAXL. .1,0.. 0..R 1ST, 0.0. -1.0.0) 

CALL  LABEL ( 1. 2+SAXL, -. 7+PAXL,. 1,0. , 0. . RLOU, 0.0. -1,0.0) 

CALL  PLOT(0. 6.0. 4,-3) 

JREC-3 

KOUNT-0 

RESET  PAXL  TO  ORIG.  VALUE  FOR  COMPUTATION  (2647  CORRECT) 

PAXL -1.25*P AXL 
SUSPEND 

CALL  EXEC(6, 0,1) 

CHECK  FOR  END  OR  PLOT  REINITIALIZATION 
IF ( ISTOP.EQ. 1 .OR. IST0P.EQ.3)  GO  TO  99 
I F ( ISTOP .EQ.-2) GO  TO  98 
RE1«R( IBUF2( 12) ) 

RE2«R(  IBUF2( 10) ) 

RE3*R(  IBUF2(  ID) 

DAT (2) -C 1*(RE2-C2) 

DAT (3) “C3+(C4*RE3) 

DAT ( 1 ) -C5*RE 1 * ( 1 + (C6*DAT ( 2 ) ) ) 

GET  TEMPERATURE  SIGN 

****  CTD  HAS  PROBLEM  UITH  SIGN  BYTE  **** 

ISGN 1 ■ I AND ( IBUF2( 13) ,2B) 

IF  ( ISGNl  .EQ  .2)  DAT(  3)—  DATC  3) 

COND-DAT1.  PRESS-DAT2,  TEMP-DAT3 

CALL  SAL  IN  CONVERSION  SUBROUTINE.  REPLACE  COND  UITH  SALIN  IN  DAT1 
CALL  LRSAL (DAT( 1 ) , DAT<2) . DAT(3) , S) 

IF ( IANS1 . EQ . 2) UR ITE ( 1.779)DAT(3) 

FORMAT ( "TEMP ■ ",F 10 .4) 

DAT  Cl) — S 

DAT1-SAL  DAT2-PRESS  DAT3-TEMP 
IF  REQUIRED,  CONVERT  TEff3  INTO  POTENTIAL  TEMP 
IF ( IANS  1 . EQ . 2) CALL  P0TMP(DAT(2) , DAT (3) , DAT( 1 ) , PTMP) 

IF ( IANS  1 . EQ . 2) DAT(3) «PTMP 
CALL  PLT(DAT ( 1 ) , DAT(3) .DAT (2) ) 

GO  TO  30 
ISTOP  — 1 
GO  TO  7 
END 

SUBROUTINE  PLT(S.T.P) 
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0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 

0111 

0112 

0113 

0114 

0115 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 


NRL  REPORT  8444 

COMMON  IREC, ISTOP, ISCAN.NSCAN, IFLOR 
COMMON  IBUFK456),  IBUF2(456) 

COMMON  C1,C2,C3.C4,C5,C6, INSTR 

COMMON  ISHIP(6), 1ST, IYR, IMO, IDAY.LOU 

COMMON  PAXL.PIV,  PINCR.TAXL.TIV, TINCR.SAXL.SIV.SINCR 

URITE(1,40)S,T, P 

URITEC6,40)S,T,P 

40  FORMAT(3F10.4) 

SMAX-S IV+CS INCR*SAXL) 

TMAX-TIV+(TINCR*TAXL) 

PMAX-P IV+CP INCR*PAXL) 
IF(S.LT.SlV.OR.S.GT.SMAX)RETURN 
IF (T. LT.TIV. OR. T.GT.TMAX) RETURN 
IF(P.LT.PIV.OR.P.GT.PMAX) RETURN 
SSCAL-SAXL/CSMAX-S IV) 

SPLT»SSCAL*(S-SIV) 

TSCAL-TAXL/CTMAX-TIV) 

TPLT-TSCAL*(T-TIV) 

PSCAL-PAXL/CPMAX-PIV) 

PPLT-PSCAL*CPMAX-P) 

NUMB1-1 

NUMB2“3 

C  ADJUST  Y  SCALE  FOR  2647  TERMINAL 

PPLT-0.8*PPLT 

CALL  SYMBISPLT.PPLT, .03.NUMB1.0.0.-1) 

CALL  SYMBCTPLT.PPLT. . 03, NUMB2, 0 . 0, -  1 ) 

RETURN 

END 

ENDS 


FTN4.L 

PROGRAM  CTDP2 
C 

C  REAL  TIME  CTD  PLOTTING  PROGRAM  UHEN  CTD  USED  IN  TOUING  MODE. 

C  PLOTS  S.T.P  VS.  TIME.  »LOT  IS  REINITIALIZED  EVERY  30  MINUTES. 

C  RELOCATABLES :  XLRSAL.  XR.  XPOTMP 

C 

C  PROGRAMMER:  LARRY  ROSENBLUM 

C  VERSION:  SEPT  18,  1979 

C 

DIMENSION  DAT(3)  .NAMES (2) , NAMETC2) ,NAMEP(2) ,NASTA(2) ,NAL0LK2) 
DIMENSION  NAMTMC4)  ,NALBL (5) , IPRAMC5) ,NAPOT (4) 

COMMON  IREC. ISTOP. ISCAN.NSCAN, IFLOR 
COMMON  IBUF 1 (456) , IBUF2 (456) 

COMMON  Cl. C2,C3,C4,C5,C6, INSTR 
COMMON  ISHIPC6). 1ST, IYR, IMO. IDAY.LOU 
COMMON  PAXL.PIV, PINCR.TAXL.TIV.TINCR.SAXL.SIV.SINCR 
DATA  NAMEP/2HPR.2HES/ 

DATA  NAMES/2HSA, 2HL  ✓ 

DATA  NAMET/2HTE,  2HMP/' 
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0021 

0022 

0023 

0024 

0025 

0026 

002? 

0028 

0029  C 

0030 

0031 

0032 

0033 

0034 

0035 

0036  7 

003?  C 
0030 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 
0058 
0059 
0060 
0061 
0062 
0063 
0064 
0065 
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
0074 
0075 
0076  C 
007?  30 

0078 
0079 


DATA  NAP0T/2HP0.2HT  -2HTE.2HMPX 
DATA  NArmV2HTI.2HrE.2H  M.2HIN/ 

DATA  NASTA/2HST.2HA  ✓ 

DATA  NAL0U/2HL0.2HU  ✓ 

DATA  NALBL/2HST.2HAR.2HT  , 2HTI.2HME/ 

CALL  RMPAR ( IPRAM) 

IANS  1  •  IPRAIi(  1 ) 

CALL  EXEC (22.8) 

FIX  AXES  LENGTH  AND  SET  TIME  AXIS  PARAMETERS 
SAXL-5. 

PAXL-5. 

TAXL-5. 

TMAXL-6.0 
TMIV-0.0 
TMINCR-5.0 
CALL  PLTLU(7) 

DRAW  AXES.  LABE,,  AND  ANNOTATE 
LR-1 

NSTIC* INT  (SAXL) +1 
NPT IC“ 1NT (PAXL) + : 

NTTlC-INTCTAXLJ+1 

KNUM1-42 

KNUM2-45 

KNUM3-46 

LNUM1 “5 

LNUM2-0 

LNUM3-3 

CALL  AX(1. 0,0. 4.6. .7, 4,0.0) 

CALL  LABEL U .0.0.2. . 1 .6. , ?, TMIV, TMINCR, -1,0.0) 
CALL  SYMB (3. 4, 0.0,  .  l.NAMTM.0.0.8) 

CALL  AX(1. 0,4. 4.4. 0,6, 4,-90. 0) 

CALL  LABEL (0.9, 4. 4, . 1,-4. ,6.PIV.PINCR.2,90. ) 

CALL  SYMB(0.9.0.0, . l.KNUMl , 90 . . - 1 ) 

CALL  SYMB (0. 1.0.5, . 1 .LNUM1 ,90 . , -  1 ) 

CALL  SYMB(0. 15.0.6. . I.NAMEP.90. ,4) 

CALL  LABEL (0 . 7,4.4, . 1 , -4. ,6, S IV, S INCR.2,90 . ) 

CALL  SYMB (0 . 7,0.8, . 1 , KNUM2, 90 . , -  1 ) 

CALL  SYMB(0. 1,2.0, . 1 , LNUM2, 98 . , -  1 ) 

CALL  SYMB(0. 15,2. 1. . I, NAMES, 90. ,4) 

CALL  LABEL (0.5, 4. 4, . 1, -4. ,6, T IV, T INCR.2,90 . ) 

CALL  SYMB (0.5, 0.0, . 1.KNUM3, 90 . ,- 1 ) 

CALL  SYMB (0. 1,3.5, . I.LNUM3.90. ,-I) 

IF ( IANS  1. EQ. 1) CALL  SYMB (0 . 15, 3 . 6, . 1 ,NAMET,90. , 4) 
IF ( IANS  1 . EQ . 2) CALL  SYMB (0 . 15,3.6, . 1.NAPOT.90.  ,8) 
CALL  SYMB (2 .5, 4. 6, . 1.NALBL.0.0, 10) 

CALL  SYM8(7. 1-3.6. . l.NASTA.0.0,4) 

CALL  SYMB(7. 1,3.3, . 1,NALOU,0.0,4) 

CALL  AX(? , 0,0 . 4, 4. , 6, 4, 90 . ) 

RIST«R( 1ST) 

RLOU'R (LOU) 

CALL  LABEL (7. 6, 3. 6, . 1 , 0. . 0 . ,R 1ST, 0.0, -  1 ,0. 0) 

CALL  LABEL (7. 6, 3. 3, . 1,0. ,0. .RLOU.0. 0,-1, 0.0) 

CALL  PLOT( 1 .0,0. 4,-3) 

JREC-3 
KOUNT  e0 
SUSPEND 

CALL  EXEC (6,0. 1) 

IF ( ISTOP.EQ. 1 .OR. I STOP. EQ. 3)  GO  TO  99 
IF ( I STOP . EQ . -2) GO  TO  98 
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0080 

0081 

0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 


IF(LR.NE. 1)  GO  TO  32 
NDAT1-IBUF2(2) 

NDAT2-IBUF2(3) 

NTMl-IBUF2(2)/24 
NTM2- IBUF2C2) -NTM1#24 
HTM3-IBUF2(3)/60 
NTM4- IBUF2 (3) -NTM3*60 
NDAT2-NDAT2-NTM4 
NTM4-0 

CALL  NUMBC3.4.4.2,  .  l,FLOATCNTMl),0.0.-l) 

CALL  NUMB(3.9.4.2, . 1, FLOAT (NTM2) ,  0.0.-1) 

CALL  NUMBC4.3.4.2, . l,FLOAT(NTM3) , 0.0.-1) 

CALL  NUM8C4.7.4.2, . I.FL0ATCNTM4) .0.0.-1) 

32  NDEL1-IBUF2(2)-NDAT1 
NDEL2-IBUF2(3)-NDAT2 
IF (NDEL l . NE.0) NDEL2-NDEL2+3600 
IF(NDEL2.GT. 1800)  GO  TO  98 
RE1-RC IBUF2C  125 ) 

RE2-RC IBUF2( 10) ) 

RE3-RC  IBUF2C  ID) 

DAT (2) -C 1#(RE2-C2) 

DAT(3)-C3+(C4*RE3) 

DAT ( 1 ) ■C5*RE 1*( 1+(C6*DAT(2) ) ) 

C  CHECK  TEMPERATURE  SIGN 

C  ISGN 1 “  IAND(  I8UF2( 12) , 28) 

C  ****  CTD  HAS  TROUBLE  IN  SIGN  BYTE  **** 

C  IF(ISGN1.EQ.2)DAT(3)—  DATC3) 

C  DAT1«C0ND,DAT2 “PRESS ,DAT3“TEMP  -  ALL  CONTVT  TO  ENG  UNITS 

C  CALL  SAL  IN  CONVERS  SUBROUTINE 

CALL  LRSAL(DAT(1),DATC2),DATC3).S) 

DAT  < 1 ) “S 

C  IF  IANS  1-2  USE  POTENTIAL  TEMPERATURE 

IF ( IANS  1 . EQ . 2) CALL  POTMP <DAT(2) , DAT(3) , DAT( 1 ) . PTMP) 

IF <  IANS1 . EQ . 2) DATC3) -PTMP 
21  CALL  PLTCDAT ( 1 ) , DAT(3) . DATC2) , NDEL2, LR) 

LR-LR+1 
GO  TO  30 

98  ISTOP— 1 
GO  TO  7 

99  END 

SUBROUTINE  PLTCS, T,P,NDEL2.LR) 

COMMON  IREC, ISTOP. ISCAN.NSCAN. IF LOR 
COMMON  IBUFK456),  IBUF2C456) 

COMMON  Cl.C2.C3.C4.C5.C6. INSTR 
COMMON  ISHIPC6), 1ST, IYR. IMO, IDAY.LOU 
COMMON  PAXL.PIV.PINCR.TAXL.TIV.TINCR.SAXL.SIV.SINCR 
UR  I TE  ( 1 , 40 )  NDEL2 ,  S ,  T,  P 
40  FORMAT( I8.3F10.4) 

SMAX-S I  V+(S  INCR>KSAXL) 

TMAX-TIV+(TINCR*TAXL) 

PMAX-P IV+(P INCR*PAXL) 

XIV-60.»0.0 

XINCR-60.*5. 

XMAX-X I V+X I NCR*6 . 

IF(S.LT.SIV.OR.S. GT. SMAX) RETURN 
IF (T.LT.TIV. OR. T.GT.TMAX) RETURN 
IF  CP . LT. P I V. OR . P . GT. PMAX) RETURN 
SSCAL-SAXLA  SMAX-S  IV) 

SPLT-SSCAL*(SMAX-S) 
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0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

014? 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

015? 

0158 

0159 

0160 


0001 

0002 

0003 

0004 

0005 

0006 

000? 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

001? 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 


TSCAL=TAXL/(TMAX-TIV) 

TPLT-TSCAL*(TMAX-T) 

PSCAL =P AXL/ ( PMAX-PIV) 
PPLT«PSCAL*(PMAX-P) 

XDEL-RCNDEL2) 

XSCAL=6 . /(XMAX-XIV) 
XPLT=XSCAL*(XDEL-XIV) 

NUMB  1=0 
NUMB2=3 
NUMB3=5 
SI2E-0. 1 

IF(LR.NE.  DSIZE-0.03 
C  SCALE  V  AXIS  FOR  264?  TERMINAL 

SPLT»0.8*SPLT 
TPLT-0.8*TPLT 
PPLT=0.8*PPLT 

CALL  SYMB ( XPLT, SPLT, SIZE, NUMB  1,0. 0,-1) 

CALL  SYMB (XPLT, TPLT, S IZE, NUMB2, 0 . 0 , -  1) 

CALL  SYMB (XPLT, PPLT,S IZE, NUMB3,0. 0,-1) 

RETURN 

END 

END* 


F-m4.L 

PROGRAM  CTDCN(3, 14) 

C 

C  CONRTOL  PROGRAM  TO  CAUSE  END  OR  IDLE. 

C  ISTOP  CHANGED  IN  COMMON  BY  *RU,CTDCN..N  WHERE 

C  N— 1  NORMAL  CONDITION 

C  N=  1  NORMAL  END 

C  i'=  2  IDLE  MODE 

C  N-  3  ABORT 

C  N— 2  BEGIN  NEU  PLOT 

C  I8UFU2(5)  USED  TO  MARK  SPECIAL  FILE  POINTS 

C  USING  IPRAM(l) 

C 

C  PROGRAMMER:  LARRY  ROSENBLUM 

C  VERSION:  SEPT  22  1979 

C 

DIMENSION  IPRAM(5) 

COMMON  IREC, ISTOP, ISCAN,NSCAN, IFLOR 
COMMON  IBUF 1 (456) , IBUF2 (456) 

COMMON  C 1 . C2.C3.C4, C5,C6, INSTR 
COMMON  ISHIP(6). 1ST, IYR, IMO, I DAY, LOU 
COMMON  PAXL,PIV,PINCR.TAXL,TIV.TINCR,SA>C-»SIV,SINCR 
CALL  RMPAR ( IPRAM) 

CALL  EXEC (22, 0) 

IF ( IPRAM(2) , EQ . 0) IPRAM<2) —1 
ISTOP- IPRAM(2) 

IBUF l (5) • IPRAM( 1) 

IBUF2(5)-IPRAM(1) 

END 

END* 
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0001 

FTN4, 

L 

0002 

PROGRAM  CTDMS 

0003 

DIMENSION  IPRAM(5) , IBACK(5) 

0004 

CALL  RMPAR (IPRAM) 

0005 

IF ( IPRAMt 1) . EQ .2)G0  TO  20 

0006 

IF ( IPRAMC 1 ) . EQ. 3)G0  TO  30 

000? 

5 

URITEC1, 10) 

0008 

10 

FORMAT ( "PLOT  TYPE?  1-S.T  VS.  P  (LOWERING) *,/ 

0009 

1  *  2-S.T.P  VS  TIME  (TOWING)*) 

0010 

READ(1.*) IANS 

0011 

IF ( IANS .NE .  1  .AND . IANS .NE . 2)  GO  TO  5 

0012 

12 

URITEd.  15) 

0013 

15 

FORMAT ( "TEMPATURE  TYPE:  1-TEMP..  2-POTENTIAL 

TEMP.  *) 

0014 

READC1.*) IANS  1 

0015 

IF ( IANS1 .NE . 1 .AND. IANS1 ,NE.2)G0  TO  12 

0016 

IBACK (2) • IANS 

001? 

IBACK (3) -  IANS  1 

0018 

COLL  PRTN( IBACK) 

0019 

CALL  EXEC(6) 

0020 

20 

WRITE( 1,25) 

0021 

25 

FORMAT ( "PAUSE  -  TYPE:  'GO.CTDRT'  FOR  DATA  AQUIS1T10N*) 

0022 

GO  TO  99 

0023 

30 

WRITE! 1,35) 

0024 

35 

FORMAT ( "PAUSE  -  IDLE  MODE  -  TYPE:  'GO.CTDRT' 

TO  CONTINUE*) 

0025 

99 

END 

0026 

END* 

0001  FTN4.L 

0002  PROGRAM  CTDFL (3, 16) 

0003  C 

0004  C  THIS  PROGRAM  SETS  THE  FLOROMETER  ON/OFF  SUICTH. 

0005  C  IFLOR-0  FLOROMETER  OFF 
0006  C  IFLOR-1  FLOROMETER  ON 

000?  C  TYPE:  'RU, CTDFL, N'  WHERE  N»2  FOR  OFF  OR  N-3  FOR  ON 
0000  C 

0009  C  PROGRAMMER :  LARRY  ROSENLUM 
0010  C  VERSION:  SEPT.  20.  19?9 
0011  DIMENSION  IPRAMC5) 

0012  COMMON  IREC, ISTOP, ISCAN.NSCAN, IFLOR 

0013  COMMON  IBUFK456),  IBUF2C456) 

0014  COMMON  Cl, C2.C3.C4.C5.C6. INSTR 

0015  COMMON  ISHIPC6) , 1ST, IYR, I MO, I DAY. LOU 

0016  COMMON  PAXL . P I V,  P INCR,  TAXL, TI V, TINCR, SAXL, S I V, S INCR 

001?  CALL  RMPAR ( IPRAM) 

0018  IFCIPRAM(l) .EQ. 2) IFLOR-0 

0019  IF(IPRAM(l) .E0.3) IFLOR-1 

0020  IF ( IPRAMC 1 ) .NE.2.AND. IPRAMd)  .NE.3)URITEC1,  10) 

0021  10  FORMAT ( "TYPE  'RU, CTDFL, N'  WHERE  N-2  (OFF)  OR  -3(0N)") 

0022  END 

0023  END* 
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0001 

0002 

0003 

0004 

0005 

0006 

000? 

0000 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

001? 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

002? 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

003? 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 


FTN4.L 

PROGRAM  CTDTR(3,98) 

C 

C  PROGRAM  TO  URITE  TERMINATION  FILE  AND  OPEN  PLOTTING  PROG.  SO 

C  THAT  IT  GETS  CLOSED  CORRECTLY.  TERM  IN  SCHEDULED  BY  CTDRT 

C  UHEN  ISTOP-  1  OR  3  -  SEE  CTDCN . 

C 

C  FORMAT  OF  B  UORD  TERMINATION  FILE: 

C  UORD  1  -1 

C  WORDS  2&3  TIME  IN  HOURS  AND  SECONDS 

C  UORD  4  ISTDP/3 

C  UORDS  5-8  0 

C 

C  PROGRAMMER:  LARRY  ROSENBLUM 

C  VERSION:  SEPT.  20,  1979 

C 

DIMENSION  IBUF4C8) , ITIME<5) .NAME  1 (3) ,NAME2(3) , IPRAMC5) 

COMMON  IREC, ISTOP, ISCAN.NSCAN, IFLOR 
COMMON  IBUF 1(456), IBUF2C456) 

COMMON  C1,C2,C3,C4,C5,C6, INSTR 

COMMON  ISHIPC6), 1ST, IYR, IMO, IDAY.LOU 

COMMON  PAXL. P I V, P INCR. TAXL, T1V, TINCR, SAXL, S IV. S INCR 

DATA  NAME 1/2HCT, 2HDP, 2H 1  /,NAME2/2HCT,2HDP,2H2  / 

C  GET  PLOT  VERSION  PARAMETER 

CALL  RMPAR(IPRAM) 

IANS-IPRAM(l) 

CALL  EXEC(22,0) 

C  URITE  EOF  AFTER  DATA 

CALL  EXEC(3, 8+1008) 

C  SCHEDULE  PLOTTING  PROGRAM 

IF C  IANS.EQ.  1)CALL  EXECU0. NAMED 
IF ( IANS.EQ.2)CALL  EXECC 10.NAME2) 

C  URITE  TERMINATION  FILE  AND  DBL  EOF 

IBUF4( 1)  —  1 
CALL  EXEC(ll.ITIME) 

IBUF4C2) «24*ITIME(5) +ITIME<4) 

IBUF4(3)-60*ITIME(3)+ITirE(2) 

IBUF4C4) • ISTOP/3 
IBUF4(5)-0 
IBUF4(6)-0 
IBUF4C7) -0 
IBUF4CB) “0 

CALL  EXEC(2,8, IBUF4.8) 

CALL  EXEC (3.8+100B) 

CALL  EXEC(3, 8+1000) 

END 

END* 
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0001 

FTN,L 

0002 

PROGRAM  F2E0F 

0003 

C 

0004 

r; 

VERSION  2/3/27 

c 

DEAN  CLAMONS 

0006 

c 

000,' 

c 

THIS  FVOCk'IM  FINDS  THE  FIRST  DO80LE  EOF  ON  A  TAPE 

0f!0-'H 

3 

HMD  POSH.  H:  Tilc  TVE  H'  lTR  THE  FIRST  OF  THE  PAIR 

I  HE  M67.IC”.'  '  •  '  TACK  b'f  T,c  RECORDS  UN  CASE  UE 

'■RE  ALPS  PDF  A’  id,.  EOF  ’  S )  AND  THEN  SPACE  FORUARD 

j  i 

ft  ft  I.  •  C  THE  NEXT  RECORD  . 

0 !  2 

0 

IF  iT  IS  AN  E  lJE  PACK  @"j  ONF.  RECOR  AND  ARE  DONE 

00  •. 

IF  Nur  AN  EUP  IJE  SPACE  FORWARD  ANOTHER  FILE 

0314 

c 

AND  REPEAT 

60  15 

c 

0016 

KITE  1  1  I!  :  ! 

03  :  2 

c 

THE  TAPE  MUST  HAVE  A  DG8S:.E  EOF  FOP  THIS 

23  i  H 

r 

PROGRAM  TO  uJORK.  IF  NOT  IT  UILl  R9N  AIJAV 

0019 

C 

0020 

LMT  '"3 

032! 

CALL  FXEC(3,200B+LMT) 

0322 

CALL  EXEC (3,200B+LMT) 

0323 

1 

CALL  EXEC(3, 1300B+LMT) 

6024 

CALL  EXEC (1, LMT, IBFF, 1) 

0025 

IF ( [EOF (LMT) . GE . 0)GO  TO  1 

0026 

CALL  EXEC (3, 200B+LMT) 

«0?7 

END 

...'"20 

ENDS 
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USER’S  GUIDE  TO  THE  REAL-TIME  CTD  PROGRAMS 
Preparation 

Prior  to  starting  the  CTD  data  collection  process  the  operator  should  ensure  that  the  system  is  set 
for  CTD  data.  The  operator  should: 

1.  Check  that  power  is  on  to  the  CTD  deck  unit  and  the  unit  is  displaying  reasonable  values. 

2.  Toggle  the  "hold"  switch  on  the  CTD  deck  unit  several  times  to  assure  that  the  sync  frame  is 
good.  The  sync  frame  indicator  should  display  either  15  or  240  (decimal). 

3.  Check  that  the  Computer  Applications  Group  CTD  Interface  has  been  switched  to  the  "ON" 
position. 

4.  If  the  computer  system  has  been  performing  other  functions,  assure  that  the  logical  units 
(LU>  assignments  are  correct.  The  CTD  programs  expect  that  LU-1  is  the  system  console. 
LU-6  is  the  line  printer,  LU-7  is  the  terminal  plotter.  U-8  is  the  magnetic  tape  unit,  and 
LU-10  is  the  CTD.  If  necessary,  logical  unit  assignments  should  be  changed  using  the  RTE 
operator  "LU"  command,  so  that  these  logical  units  will  point  to  the  appropriate  equipment 
table  entry. 

5.  Mount  a  tape  on  the  tape  drive  and  check  that  the  subchannel  indicator  on  the  tape  unit  is 
set  to  0  and  that  the  tape  is  "write"  enabled  (the  write  ring  is  installed). 

6.  Obtain  the  latitude  and  longitude  from  a  laboratory  display,  a  computer  program  for  this  pur¬ 
pose,  or  the  navigation  laboratory. 

If  a  CTD  lowering  (or  tow)  has  just  been  completed,  then,  to  begin  the  next  lowering  lor  tow), 
the  operator  should  get  an  RTF  system  prompt  *  and  type  "RU,  F2EOF."  This  program  will  position 
the  tape  between  the  two  HOI 's  which  follow  the  termination  record  of  the  previous  CTD.  Perform 
step  6  above. 

Program  Start-Up 

The  next  two  sections  of  this  chapter  describe  the  interaction  between  the  operator  and  CTD 
real-time  programs.  For  purposes  of  identification  within  this  report  only,  all  user  responses  are  under¬ 
lined  All  user  entries  are  to  be  followed  by  a  carriage  return.  The  symbol  *  is  the  RTE  system  prompt 
and  may  be  obtained  in  an  RTE  system  by  striking  any  key  on  the  console  one  time.  Answers  to  ques¬ 
tions  are  by  way  of  example. 

*  RU,  CTDRT 

ENTER  CTD  MODE.  I  =  LOWERING.  2  =  TOWING 

j 

FNTER  TEMPERATURE  TYPE,  1  =  TEMP,  2  =  POTENTIAL  TEMP 

2 

ENTER  SHIP’S  NAME  (CRUISE)  (UP  TO  1 2  CHARACTERS) 

HAYES  9999 

I  N  K  R  DA  II.  (YR.  MO,  DAY) 

1474  10  8 

I  NTER  STATION  NUMBER  &  LOWERING  NUMBER 

2J 

ENTER  LATITUDE:  DEG.  MIN.  &  FRACTION  OF  MIN  xxxx..  xx  xxx 
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32.  15.25 

ENTER  LONGITUDE  DEG,  MIN.  &  FRACTION  OF  MIN  xxxx.,  xx.xxx 
-65.  22.23 

At  this  point  the  operator  is  asked 

USE  DEFAULT  VALUES,  1  =  YES,  2  *=  NO 

The  default  values  concern  the  instrument  number,  the  values  of  the  constants  Cl.  C2.  ....  C 6 
used  to  convert  readings  into  pressure,  temperature,  and  conductivity  and  the  axes  limits  for  the  real¬ 
time  plots.  The  conversion  constant  default  values  are:  Cl  =  0.1.  C 2  =  0.0,  C3  =  0  0,  C4  =  0.0005, 
C5  =  0.001,  and  C6  =  0.0.  The  only  constant  with  which  the  operator  need  be  concerned  is  Cl. 
which  is  dependent  on  the  depth  rating  of  the  pressure  sensor.  The  default  is  for  the  6.500-decibar 
(65-MPa)  pressure  sensor.  If  a  different  pressure  sensor  is  used,  Cl  should  be  changed  to  the  value 
shown  in  Table  1. 


Table  I  —Values  of  the 
Conversion  Constant  C’l  for  Each  of  the 
Different  CTD  Pressure  Sensors 


Pressure-Sensor  Range 
(decibar)  (MPa) 

Cl  Value 

0-6500 

0-65 

0.1  (default) 

0-3200 

0-32 

0.05 

0-1600 

0-16 

0.025 

0-650 

0-6.5 

0.01 

The  axes  lengths  are  fixed  at  6  inches  (15.2  cm)  for  the  a  axis  and  5  inches  (12.7  cm)  for  the  .i 
axis.  The  operator  supplies  the  initial  value  and  the  increments  per  inch  for  each  parameter— pressure, 
temperature,  and  salinity.  The  default  values  for  lowering  plots  are  given  in  Table  2  and  for  towing 
plots  in  Table  3.  In  general,  except  for  deep-water  lowerings,  the  operator  will  want  to  specify  more 
limited  axes  scaling  using  knowledge  of  the  area  and  depth  involved.  Returning  to  the  question  about 
default  values,  let  us  assume  we  wish  to  lower  a  1600-decibar  (16-MPa)  pressure  sensor  from  200  to 
500  decibars  (2  to  5  MPa)  and  that  we  know  that  temperature  will  be  in  the  range  14  to  19°C  and  salin¬ 
ity  in  the  range  34.5  to  36.0  mmho(mS).  In  the  following  the  symbol  "n  CR”  will  denote  a  single  space 
followed  by  a  carriage  return.  Then,  continuing  our  sample  question  and  answer,  the  programs  would 
ask: 

USE  DEFAULT  VALUES,  1  =  YES,  2  =  NO 
2 

DEFAULT  VALUES  GIVEN  IN  PARENTHESIS 
INSTRUMENT  #  (1) 

□  CR 

CAL  CONST  PRESSURE  Cl,  C2  (0.1, 0.0) 

0.025  0.0 

CAL  CONST  TEMP  C'3,  C4  (0.0.  0.0005) 

□  CR 

CAL  CONST  C5,  C6  (0.001,  0.0) 

□  CR 

FOR  EACH  AXIS  THE  AXIS  LENGTH  IS  GIVEN  IN  PARENTHESIS  ENTER 
THE  INITIAL  VALUE  AND  THE  INCREMENT/INCH. 

PRESSURE  (AX  LNTH  ==  5)  ENTER  INIT  VALUE,  INCR  INCH 
200.0  60.0 

TEMP  (AX  LNTH  =  6)  ENTER  INIT  VALUE,  INCR/INCH 
14.0  1.0 

SAL  (AX  LNTH  =  6)  ENTER  INIT  VALUE,  INCR/INCH 
34,5  0.25 
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Table  2  — Plot  Default  Values  for  the  Lowering  Mode 


Parameter 

Axis  Length 
(in.)  (cm) 

Initial  Value 

Incr 

(per  in.)  (per  cm) 

Range 

Pressure 

5.0 

12.7 

0 

1000  394 

0-5000  dbar 
(0-50  MPa) 

Temperature 

6.0 

15.2 

1 

4  1.6 

1-25  °C 

Salinity 

6.0 

15.2 

32 

1  0.4 

32-38  mS 

Table  3  -  Plot  Default  Values  for  the  Towing  Mode 


Parameter 

Axis 

(in.) 

Length 

cm 

Initial  Value 

(per  in. 

Incr 

)  (per  cm) 

Range 

Pressure 

5.0 

12.7 

0 

10 

3.9 

0-50  dbar 

Temperature 

5.0 

12.7 

1 

4 

1,6 

(0  -  0.5  MPa) 

Salinity 

5.0 

12.7 

33 

1 

0.4 

33  -  38  mS 

Time 

6.0 

15.2 

Start  time  of  plot 

5 

2.0 

0—30  min 

The  operator  entries  are  printed  out  on  the  terminal  and  written  to  the  line  printer.  The  operator 
is  then  given  the  option  of  restarting  the  question  ptueess  if  any  mistakes  have  been  made  (the  answer 
"3"  to  the  following  question  is  for  debug  purposes  only  and  should  not  be  used). 

ARE  INPUT  PARAMS  O.K.:  1  =  YES,  2  =  NO,  3  =  YES  -  NO  MT  FILE 
1 

PAUSE-TYPE  ‘GO,  CTDRT’  FOR  DATA  ACQUISITION 
At  this  point  the  plot  axes  and  annotation  will  be  output  on  the  plotting  terminal.  The  operator  should 
verify  that  these  are  correct,  especially  the  parameter  limits  on  each  axis.  Once  the  plot  axes  and  anno¬ 
tation  are  drawn  the  program  is  prepared  to  acquire  data.  The  operator  now  obtains  the  system  prompt 
(*)  and  types: 

*  GO,  CTDRT 

This  completes  the  start-up  procedure. 

Special  Functions 

The  program  CTDCN  is  used  to 

•  end  data  taking, 

•  abort  data  taking, 

•  put  the  program  into  an  "idle"  mode,  or 

•  mark  special  segments. 


The  end  and  abort  produce  exactly  the  same  result,  except  that  a  word  in  the  termination  file  indi¬ 
cates  which  way  the  program  was  terminated.  The  idle  mode  is  used  to  suspend  data  collection  tem¬ 
porarily  and  the  special  segment  marker  changes  one  word  in  the  data  buffers.  This  last  operation  can 
be  used  to  identify  data  segments  for  later  processing.  These  operations  occur  when  the  operator  types: 
for  a  NORMAL  TERMINATION 

*  RU,  CTDCN,,  1 
for  an  ABORT 

*  RU,  CTDCN,  ,3 
for  an  IDLE 

*  RU,  CTDCN,  ,2 
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The  message  PAUSE-IDLE  MODE-TYPE:  'GO.  CTDRT’  TO  CONTINUE  will  appear  on  the  sys¬ 
tem  console.  To  continue  the  collection  process  the  operator  types: 

*  GO.  CTDRT 

for  marking  SPECIAL  SEGMENTS 

*  RU,  CTDCN.N 

where  N  is  any  number.  This  number  is  placed  in  Word  5  of  each  data  buffer  until  it  is  changed  or  the 
programs  are  terminated. 

The  program  CTDFL  sets  the  fluorometer  on/off  switch. 

The  command 

*  RU.CTDFL.  3 


turns  the  switch  on  (fluorometer  data  being  taken)  and  the  command 
*  RU.  CTDFL.  2 

turns  the  switch  off.  It  should  be  noted  that  the  switch  is  initialized  as  "off'  and  that  it  is  used  only  as  a 
marker;  failure  to  set  the  switch  does  not  alter  the  data. 

DATA-PROCESSINo  PROGRAMS 
Parameters  of  Interest 

Once  a  CTD  data  tape  has  been  collected,  many  parameters  of  interest  to  the  oceanographic 
researcher  can  be  computed  from  the  stored  values  of  conductivity,  temperature,  and  pressure.  To 
assist  the  researcher  a  program  called  PLTBK  has  been  written  which  will  compute  values  and  produce 
plots  of  various  parameters.  To  do  this  the  user  is  offered  the  following  parameter  list: 

1  =  Salinity 

2  —  Pressure 

3  —  Temperature 

4  =  Potential  Temperature 

5  =  Sound  Velocity 

6  =  Density 

7  =  Specific  Volume 

8  =  Adiabatic  Temperature  Gradient 

9  =  Sigma  (0) 

10  =  Depth 

1 1  =  Brunt-Vaisala  Frequency 

12  =  Time  (minutes) 

13  =  d(temp)  /  d( pressure) 

14  =  d(salinity)  /  d(pressure) 

15  =  Conductivity 

16  =  Fluorometer 


Program  PLTBK  allows  the  user  to  compute  and  plot  (against  each  other)  any  pair  of  the  above 
parameters  or.  alternately,  to  define  any  one  parameter  as  the  y  axis  and  plot  any  two  variables  as  .v-axis 
parameters  against  the  selected  v-axis  parameters.  Thus,  for  example,  if  the  user  wanted  to  obtain  a 
sound-speed  profile  the  user  would  select  5(sound  velocity)  for  the  v  axis  and  lO(depth)  for  the  y  axis. 
To  obtain  a  plot  of  temperature  and  conductivity  vs  pressure,  the  user  would  select  3(temperature)  for 
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the  v  axis,  2(pressure)  for  the  y  axis  and  15 (conductivity)  as  a  second  \  axis.  The  axes  scales  for  the 
plot  are  user  selected.  Scientific  subroutines  for  calculating  these  parameters  are  listed  in  the  appendix. 
These  subroutines  were  provided  by  the  Woods  Mole  Oceanographic  Institution,  except  for  the  Bruni- 
Vaisala  frequency  calculation,  which  uses  an  approximation,  formula  and  is  valid  for  the  top  several 
hundred  meters. 

User’s  Guide  to  Program  PLTBK 

This  section  describes  the  interaction  between  the  user  and  program  PLTBK.  As  before,  user 
responses  will  be  underlined  and  answers  are  by  way  of  example. 

*  RU,  PLTBK 

ENTER  STATION  NUMBER  AND  LOWERING  NUMBER 
2J 

DO  YOU  WANT  A  TIME  SEARCH ?  1  =  YES,  2  =  NO 
A  time  search  allows  the  user  to  look  at  only  part  of  an  event;  otherwise  the  entire  event  will  be  pro¬ 
cessed. 

\ 

ENTER  START  TIME-DAY  HR  MIN  SEC 
171  5  27  0 

ENTER  STOP  TIME-DAY  HR  MIN  SEC 
171  6  31  0 

1  “SALINITY  2  =  PRESSURE  3=TEMPERATURE 

4  =  POT.  TEMP  5=SOUNDVEL  6  =  DENSITY 

7  =  SPEC  VOL  8  =  ATG  9  =  SIGMA 

10  =  DEPTH  1 1  =B-V  FREQ  12  =  TIME 

13  =  D(T)/D(P)  14  =  D(S)/D(P)  15  “CONDUCTIVITY 

16=FLUOROMETER 

ENTER  X  AXIS  PARAMETER  NUMBER 
3 

ENTER  Y  AXIS  PARAMETER  NUMBER 

JO 

DO  YOU  WANT  A  SECOND  X  AXIS,  1  =YES,  2  =  NO 
\ 

ENTER  X2  AXIS  PARAMETER  NUMBER 
15 

ENTER  X  AXIS  LENGTH,  INIT  VALUE,  INCR/INCH 
6  3.  15. 

ENTER  Y  AXIS  LENGTH,  INIT  VALUE,  INCR/INCH 

5  2000.  400. 

ENTER  X2  AXIS  LENGTH,  INIT  VALUE,  INCR/INCH 

6  32.  .5 

Since  using  all  points  usually  provides  too  dense  a  plot,  the  user  chooses  the  plotting  density  by 
answering  the  following  question: 

#  OF  POINTS/4  SEC  SCAN;  MUST  DIVIDE  64 

8 

The  program  will  now  plot  the  data  and  stop. 

Examples  of  Program  PLTBK  Graphics 

In  this  section,  several  examples  of  program  PLTBK  plots  are  given.  The  plots  were  produced 
using  PLTBK  from  a  digital  tape  recorded  on  the  USNS  Lynch  in  April  1978  in  the  Atlantic  off  the  con¬ 
tinental  sheif  of  the  United  States.  Figure  6  shows  salinity  and  temperature  plotted  against  pressure  for 
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the  full  lowering  (approx.  100  m  to  4600  m)  plotting  one  point  every  4  seconds,  f  igure  7  shows  salin¬ 
ity  for  the  top  1000  meters,  plotting  one  point  every  second.  Figure  8  is  the  sound-velocity  profile  for 
the  time  and  site  Figure  9  is  a  temperature-salinity  plot  for  the  lowering.  Figure  10  is  a  density  plot 
and  Figure  11  plots  density  vs  potential  temperature  Illustrating  the  ability  to  plot  any  two  parameters 
against  any  other.  Figure  12  shows  salinity  and  temperature  plotted  against  sound  speed. 
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Figure  6  —  Plot  of  salinity  and  temperature 
vs  pressure  for  the  entire  lowering 
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Data  Processing  Programs— Fortran  Source  Listings 


0001  FTN4,  L 

0002  PROGRAM  PLTBK 

0003  C 

0004  C  THIS  VERSION  REVISED  TO  WORK  FOR  FLUOROMETER  DATA  SET  10  /73 

0005  C  THIS  PROGRAM  (JILL  PRODUCE  PLOTS  FROM  THE  COLLECTED  CTD  DATA 

0006  C  TAPE.  ANY  PARAMETER  FROM  THE  PARAMETER  LIST  CAN  BE  PLOTTED 

000?  C  AGAINST  ANY  OTHER  PARAMETER  OR  AGAINST  ANY  PAIR  OF  PARAMETERS 

0000  C  FROM  THE  LIST.  TAPE  MUST  BE  POSITIONED  AT  THE  DATA.  SCIENTIFIC 

0009  C  SUBROUTINES  USED  WERE  SUPPLIED  BY  (J.H.O.I.  EXCEPT  FOR  THE  APPROX. 

0010  C  USED  FOR  BRUNT- VA 1SALA  FREQUENCY. 

0011  C 

0012  C  RELOCATABLE  SUBROUTINES:  *LRSAL,«R,KP0TMP,*SVEL,8ATG,«SIGMT, 

0013  C  2SPV0L.KBRVA1 

0014  C 

0015  C  THIS  VERSION  ADDS  THE  FLUOROMETER. 

0016  C 

0017  C  PROGRAMMER:  LARRY  ROSENBLUM  VERSION:  SEPT.  1979 

0018  C 

0019  DIMENSION  DATC 16) .NAMEP (8) .NAMES (8) . NAMETC8) . NAMPTC8) .NAME (8. 16) 

0020  DIMENSION  NAMEV(8) ,NASTA(2) ,NAL0U(2) ,NADEN(8) ,NASPV(8) ,NAATG(8) 

0021  DIMENSION  NASMT(B) .NADPH (8) .NABVFC8) , NATME (8) ,NALBL (5) ,NADLTC8) 

0022  DIMENSION  NADLS(B) ,NAFLR(8) ,NACON(B) 

0023  DIMENSION  IREG (2) . XBVC500) , IAX13) 

0024  DOUBLE  PRECISION  XALPH , XB V, DENS , SUMRH .XSO.XST 

0025  COMMON  IREC, ISTOP, ISCAN.NSCAN. IFLOR 

0026  COMMON  I8UF 1 (456) , IBUF2 (456) 

002?  COMMON  C1.C2,C3,C4.C5,C6, INSTR 

0028  COMMON  ISHIPC6) , 1ST, IYR, IMO, IDAY,LOU 
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0029  COMMON  X1AXL.X1 IV,X1 INCR,YAXL,YIV,YINCR,X2AXL,X2IV.X2INCR 

0030  EQUIVALENCE  (NAMES, NAME ( 1, 1) ) , (NAMEP.NAME ( 1,2) ) 

0031  EQUIVALENCE  (NAMET, NAME ( 1 . 3) ) , (NAMPT. NAME( 1 ,4) ) 

0032  EQUIVALENCE  ( NAME V, NAME (1,5) ) , (NADEN.NAME( 1.6) ) 

0033  EQUIVALENCE  (NASPV, NAME (1,7)), (NAATG.NAMEC 1,8) ) 

0034  EQUIVALENCE  (NASMT, NAME (1,9)), (NADPH,NAME( 1, 10) ) 

0035  EQUIVALENCE  (NABVF,NAME( 1, ID), (NATME.NAME (1 , 12) ) 

0036  EQUIVALENCE  (NADLT, NAME (1,13)), (NADLS.NAME (1, 14) ) 

0037  EQUIVALENCE  (NACQN,NAME( 1, 15) ) , (NAFLR,NAME( 1, 16) ) 

0038  EQUIVALENCE  (REG, IREG) 

0039  DATA  NAMEP/2HPR, 2HES, 2HSU,2HRE, 2H  ,2HDB,2HAR,2HS  / 

0040  DATA  NAMES/2HSA , 2HL 1 , 2HN 1 , 2HTY, 2H  ,2H  .2H0/.2H00/ 

0041  DATA  NAMET/2HTE.2HMP, 2HER.2HAT, 2HUR.2HE  , 2HDG.2H  C/ 

0042  DATA  NAMPT/2HP0, 2HTE, 2HNT.2H IA, 2HL  ,2H  ,2HTE,2HMP/ 

0043  DATA  NAMEV/2HVE , 2HLO, 2HC I , 2HTY, 2H  ,2HM/,2HSE,2HC  / 

0044  DATA  NADEN/2HDE , 2HNS, 2H I T, 2HY  , 2H  G, 2HM/, 2HCM, 2H3  / 

0045  DATA  NASPV/2HSP , 2H .  ,2HVO,2HL. ,2H  ,2HCM.2H3/,2HGM/ 

0046  DATA  NAATG/2HAD, 2H IA, 2HBT, 2H  T.2HEM.2HP  , 2HGR, 2HAD/ 

0047  DATA  NASHT/2HSI.2HGM.2HA  ,2H  T.2H  ,2HGM,2H/C,2HM3/ 

0040  DATA  NADPH/2H  ,2HDE.2HPT,2HH  ,2H  M, 2HET, 2HER, 2HS  / 

0049  DATA  NABVF/2H8R,2H-V,2HA  ,2HFR,2HEQ,2H  C.2HY/.2HHR/ 

0050  DATA  NATME/2H  , 2HTI. 2HME, 2H  ,2H  .2HMI.2HN  ,2H  / 

0051  DATA  NADLT/2H  ,2HDE,2HL(,2HT),2H/D,2HEL,2H(P,2H)  / 

0052  DATA  NADLS/2H  ,2HDE.2HL(,2HS),2H/D,2HEL,2H(P,2H)  / 

0053  DATA  NAC0N/2H  ,2HC0.2HND.2HUC,2HTI,2HVI,2HTY.2H  / 

0054  DATA  NAFLR/2H  ,2H  F, 2HL0, 2HR0, 2HME, 2HTE, 2HR  ,2H  / 

0055  DATA  NASTA/2HST.2HA  / 

0056  DATA  NAL0U/2HL0, 2HU  / 

0057  DATA  NAL8L/2HST, 2HAR, 2HT  .2HTI.2HME/ 

0058  C 

0059  LUMT-8 

0060  POLD-0.0 

0061  URITE(l.lll) 

0062  111  FORMAT  ( •■ENTER  STATION  NUMBER  AND  LQUERING  NUMBER") 

0063  READ(1,*) 1ST. LOU 

0064  IJR  ITE(  1,200) 

0065  200  FORMAT! "DO  YOU  UANT  A  TIME  SEARCH?  1-YES.  2-NO") 

0066  READ! 1,*) I ANSI 

0067  IF(IANSl.NE.l)  GO  TO  206 

0068  URITE( 1.202) 

0069  202  FORMAT! “ENTER  START  TIME  -  DAY  HRS  MIN  SECS") 

0070  READ (1, >10  IST1,  IST2,  IST3,  IST4 

0071  UR ITE( 1,204) 

0072  204  FORMAT( "ENTER  STOP  TIME  -  DAY  HRS  MIN  SECS") 

0073  READ(1.*) ISP1, ISP2, ISP3, ISP4 

0074  !SP5=24*ISP1+ISP2 

0075  IST5-24*ISTl+IST2 

0076  IST6-60*IST3+IST4 

0077  ISP6-60*ISP3+ISP4 

0078  206  UR ITE( 1 , 12) 

0079  12  FORMATC  1-SALINITY  2-PRESSURE  3-TEMPERATURE",/' 

0080  I  "  4-POT.  TEMP  5  -SOUND  VEL.  6-DENSITY",/ 

0081  2  "  7-SPEC.  VOL.  8-ATG  9-SIGMA  T(MOD) ",/ 

0082  3  “10-DEPTH  Il-B-VFREQ.  12-TIME(MIN) "./ 

0083  4  “13*D(T)/D(P)  14-D(S)/D(P)  15-CONDUCTIVITY"./ 

0084  5  " 16-FL0R0METER  *  > 

0085  URITE( 1, 101) 

0086  101  FORMAT( "ENTER  X  AXIS  PARAMETER  NUMBER") 

0087  READ(1.*)  IAXd) 
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0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

009? 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

010? 

0108 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

011? 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

012? 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

013? 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 


URITEd,  102) 

102  FORMAT! “ENTER  Y  AXIS  PARAMETER  NUMBER") 

READ! 1,*) !AX(2) 

URITE! 1, 105) 

105  FORMAT! “DO  YOU  UANT  A  2ND  X  AXIS:  1=YES.  2=N0") 

READ! 1,*) I OPT 

IF! I0PT.EQ.2)  GO  TO  123 
URITE! 1, 106) 

106  FORMAT! “ENTER  X2  AXIS  PARAMETER  NUMBER") 

READ! 1, HO  IAX!3) 

123  IF!IAX(l).NE.ll .AND. I AX! 3) -HE. 1 1 .AND. IAX!2) .NE. 1 1 )GO  TO  113 
URITEd, 12?) 

12?  FORMAT! “BE  SURE  THAT  B-V  FREQ.  IS  PLOTTED  VS.  DEPTH  ONLY") 
URITE! 1, 124) 

124  FORMAT !  “ENTER  B-V  FREQ.  DEPTH  BIN  SIZE,  BETUEEN  111.  8,  20  M.“) 
READ!  l,HOAVDPT 

IF ! AVDPT. LT. 1 . 0 . OR . AVDPT. GT. 20 . 0) UR ITE ! 1 , 1 26) 

126  FORMAT! “FREQ.  BIN  NOT  UITHIN  REQUIRED  RANGE") 

IF ! AVDPT.LT. 1 . 0. OR. AVDPT. GT. 20.0) GO  TO  123 
URITEd,  125) 

125  FORMAT! "ENTER  SURFACE  DENS1TY-OBTAIN  FROM  DENSITY  PLOT") 

READ!  I.XORHO0 

C  INITIALIZE  BRUNT- VA IS ALA  AVERAGING  COUNTERS 

NUMBV- 1 

DPTBV=0 . 0+AVDPT 
DLAST=0 . 0 

113  IF! IAX! 1 ) .NE. 12. AND . IAX!3) .NE. 12. AND . I AX! 2) .NE . 12 . )  GOTO  10? 
URITE! 1, 110) 

110  FORMAT! “ENTER  PLOT  START  TIME:  DAYS, HRS, MIN, SECS") 

READ! 1,*)NT1,NT2,NT3,NT4 
NDAT1=24*NT1+NT2 
NDAT2=60*NT3+NT4 
10?  UR ITE! 1,103) 

103  FORMAT! "ENTER :X  AXIS  LENGTH,  IN  IT.  VALUE.  INCR/INCH") 

READ ! 1 , *) XI AXL. XI IV, XI INCR 

URITE! 1, 104) 

104  FORMAT! “ENTER:  Y  AXIS  LENGTH,  INIT.  VALUE,  INCR/INCH") 

READ! 1,*) YAXL, YIV, YINCR 

IF! IOPT.EQ.2)  GO  TO  109 
URITE! 1, 108) 

108  FORMAT! “ENTER:  X2  AXIS  LENGTH,  INIT.  VALUE.  INCR/INCH") 

READ! 1, *) X2 AXL , X2 I V, X2 I NCR 

109  CALL  IN IT1 
URITE! 1, 11) 

11  FORMAT!"*  OF  PTS/4  SEC  SCAN;  MUST  DIVIDE  2**?“) 

READ! 1,*) INUM 
JSK IP =448/1 NUM 
CALL  PLTLU! ?) 

NYTIC=  INT(YAXl.)  +1 

NX1TIC=INT!X1AXL)+1 

NX2TIC=INT!X2AXL)+1 

CALL  AX!0.6,0.4,XIAXL,NXITIC,3,0.0) 

CALL  AX!0 . 6, YAXLH0 . 4, YAXL.NYTIC, 3, -90 . ) 

CALL  LABEL(0.6,0.2. . 1,X1AXL,NX1TIC.X1 IV.Xl INCR, 3, 0.0) 

CALL  SYMB1-.3+X1AXL/2. ,0.0, . 1 , NAME ( 1 , IAX! 1)), 0.0, 16) 

CALL  LABEL !0 . 3, YAXL+0 . 4, . 1 ,-YAXL, NYTIC . YIV, YIHCR, 2, h  90 . ) 

CALL  SYM8 !0 . 1 , - . 3+YAXL/2 . . . 1 , NAME ! 1 , IAX!2) ) ,90.0, 16) 

CALL  S YMB ! . 7  +X 1 AXL , - . 4+YAXL , . 1,NASTA,0.0,4) 
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0146 
014? 
0148 
0149 
0150 
0151 
0152 
0153 
0154 
0155 
0156 
015? 
0158 
0159 
0160 
0161 
0162 
0163 
0164 
0165 
0166 
016? 
0168 
0169 
01?0 
0 1?  1 
0 1?2 
0 !  ?3 
0 1?4 
0175 
0176 
017? 
0178 
0179 
0180 
0181 
0182 
0183 
0184 
ei85 
0186 
018? 
0188 
0189 
0190 
0191 
0192 
0193 
0194 
0195 
0196 
0197 
0198 
0199 
0200 
0201 
0202 
0203 
0204 


CALL  S  YMB ( . 7+X 1 AXL , - . 7 +YAXL , . l.NALCU,0.0,4) 

RIST-RC 1ST) 

RLOU-RCLOU) 

CALL  LABEL! 1 . 2+X1AXL.-.4+YAXL, . 1, 0. , 0. ,R IST.0 .0,- 1 .0.0) 

CALL  LABEL (1.2+XIAXL.-.7+YAXL. . 1,0. ,0. ,RLOU,0. 0.-1, 0.0) 

IF ( 10PT.EQ.2)  GO  TO  19 

CALL  AXC .6+X2AXL, .4+YAXL,X2AXL,NX2TIC,3,-180. ) 

C ALL  L ABEL ( . 6 , . 6  +YAXL . . 1 , X2AXL , NX2T I C , X2 1 V, X2 1 NCR , 2, 0 . ) 

CALL  SYM8C-.3+X2AXL/2., .8+YAXL, . I , NAME ( 1 , IAXC3) ) ,0. 0, 16) 

19  CALL  PLOT(0. 6,0. 4,-3) 

IF ( IAXC l ) .NE . 12 . AND . IAXC2) .NE . 12. AND. IAX(3) . NE . 12)  GO  TO  15 
CALL  SYMBC2. 1,5.2, . 1,NALBL,0.0, 10) 

CALL  NUMBC3.4.5.2, . 1 ,FLOATCNTl ) ,0.0, -1 ) 

CALL  NUhB(3.9,5.2, . 1 , FLOAT (NT2) ,0. 0, -1 ) 

CALL  NUMB (4. 3, 5. 2, . 1, FLOAT (NT3) .0.0.-1) 

CALL  NUM8(4.7,5.2, . 1, FLOAT (NT4) ,0.0,-l) 

C  READ  IN  DATA 

15  REG-  EXECC1, LUMT, IBUF2, 456) 

ILEN-IREGC2) 

IF ( IANS  1 .NE . 1)  GO  TO  208 
IF ( IEOF (LUMT) .LT.0)  GO  TO  15 
IF ( IBUF2(2) .LT. IST5)  GO  TO  15 
IF( ILEN.NE.456)  GO  TO  15 

IF( IBUF2(2) .EQ. IST5.AND. IBUF2C3) .LT. ISTG)  GO  TO  15 
IF ( IBUF2(2) .GT. ISPS) GO  TO  99 

IFC IBUF2C2) .EQ. ISP5. AND. IBUF2C3) .GT. ISP6)  GO  TO  99 
208  IF ( IEOF (LUMT) .LT.0)  GO  TO  99 
DO  20  1-9,456, JSK IP 
RE  1 =R ( IBUF2 ( 1+3) ) 

RE2-RC IBUF2( 1+1) ) 

RE3-RC IBUF2( 1+2) ) 

C 

C  CALCULATE  CONDUCTIVITY,  PRESSURE.  AND  TEMPERATURE 

C  DAT1-COND.  DAT2-PRESS,  DAT3-TEMP 

C  THEN  CONVERT  TO  SALINITY  AND  PUT  DAT 1 -SAL 

C  STORE  COND  IN  DAT(15) 

C 

DAT(2)“Cl*(RE2-C2) 

DATC3)=C3+CC4*RE3) 

DAT(  1)  -C5*RE  1*(  1+(C6*DAT(2) ) ) 

C  GET  SIGN  OF  TEMPERATURE 

C  ****  MAY  BE  PROBLEM  IN  CTD  WITH  STATUS  BIT  **** 

C  ITST1-IAND( IBUF2( 1+4), 2B) 

C  IFC ITST1.EQ.2)  DAT(3) — DATC3) 

DAT( 15) -DAT (l) 

CALL  LRSAL  (DA1'(  1) ,  DATC2)  ,DAT(3)  ,S) 

DAT(1)=S 

URITE(6,992)DAT(2) ,DAT(3) ,DAT( 1) 

992  FORMHT( IX, 3CF12.2, 2X) ) 

C 

C  CALCULATE  FLOROMETER  READING 

IFC IAX(1) .NE. 16. AND. IAXC2) .NE. 16. AND. IAXC3) .NE. 16)G0  TO  21 
DAT ( 16) -FLOAT ( 1BUF2C 1+5) ) *FLOAT( IBUF2C 1+6) ) 

C 

C  CALCULATE  POTENTIAL  TEMPERATURE 

C 

21  IF C IAXC 1 ) . NE . 4. AND . IAXC 2) .NE . 4. AND. IAXC3) .NE.4)  GO  TO  24 
CALL  P0TMPCDATC2) ,  DAT! 3) , DAT ( 1 ) ,PTMP) 

DATC4) -PTMP 
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0205  C 

0206  C  CALCULATE  SOUND  VELOCITY 

020?  C 

0209  24  IF  t IAX( 1) .NE.5.AND. IAX(2) .NE.5.AHD. IAX(3) .NE.5)  GO  70  26 

0209  CALL  SVEL (DAT (2) , DAT (3) , DAT( 1 ) , DAT (5) ) 

0210  C 

0211  C  CALCULATE  SPECIFIC  VOLUME  AND  DENSITY 

0212  C 

0213  26  IF ( IAX(l) .EQ.7.0R. IAX(3) .EQ.7.0R. IAX(2) .E0.7)  GO  70  29 

0214  IF(IAX(1)  .EQ.6.0R.  IAXC3)  .EQ.6.0R.  IAX(2)  .EQ.6)  GO  TO  29 

0215  IF(IAXd)  .EQ.9.0R.  IAX(3)  .EQ.9.0R.  IAX(2)  .E0.9)  GO  TO  29 

0216  I F C I AX C 1 1 .EQ . 1 1 . OR . IAX(3) . EO. 1 1 . DR . I AX(2) . EU . 1 1 )  GO  U  29 

021?  GO  TO  28 

0218  29  CALL  S IGTIT ( DOT C3 ) ,  DATC  1) ,  XS0,XST) 

0219  CALL  SP VOL (DAT (2),DAT(3). XS0,XST, XALPH) 

0220  DAT (?) = XALPH 

0221  DENS=1.0/XALPH 

0222  DAT (6) = 1 . /DAT (7) 

0223  DAT (9) *  1000 .*(DAT(6)-1 .0) 

0224  C 

0225  C  CALCULATE  ADIABATIC  TEMPERATURE  GRADIENT 

0226  C 

022?  28  IFCIAX(l) .NE.8.AND. 1AX(3) .NE.8.AHD. IAXC2) .NE.9)  GO  TO  30 

0228  CALL  ATG(DAT(2) , DAT(3)  ,  DAT <  1 ) , GAMMA) 

0229  DAT(8) “GAMMA 

0230  C 

0231  C  CALCULATE  DEPTH  FROM  PRESSURE 

0232  C 

0233  30  IF ( IAX( 1 ) .EQ . 10.OR . IAX(2) .EQ. 10 .OR . IAX(3) .EQ . 10)  GO  TO  32 

0234  IF ( IAX( 1 ) . EQ . 1 1 . OR . IAXC2) .EQ. 1 1 .OR . IAX(3) .EQ . 1 1)G0  TO  32 

0235  GO  TO  31 

0236  32  IF (DAT (2) . LE. 600.0) DAT ( 10) “DAT(2)/1 .009 

0237  IF(DAT(2) .LE.600.0)GO  TO  31 

0238  A=. 980665*2. 022365 

0239  A2-1.E-6 

0240  A-A*A2 

0241  8 *. 980665*1. 02940 

0242  C--DAT(2) 

0243  ROOTZ=>SQRT ( (8**2)  -4*A*C) 

0244  DAT ( 10) * (-B+ROOTZ) /(2*A) 

0245  C  PRES(DBARS) “.980665*( 1 .02940Z+2.022365*10**-6*Z**2)  Z  IN  METERS 

0246  C 

0247  C  CALCULATE  BRUHT-VAISALA  FREQUENCY 

0248  C  AVERAGE  DENSITIES  WITHIN  EACH  BIN  AND  THEN  COMPUTE  BV 

0249  C 

0250  31  IFCIAXCI) .NE. 11 .AND. IAX(2) -NE. 11 .AND. IAXC3) .NE. 1 1)  GO  TO  756 

0251  C  CHECK  FOR  1ST  POINT  OF  NEW  SEGMENT 

0252  IFCNUMBV.NE. 1)  GOTO  33 

0253  DMXZ'DPTBV 

0254  DMN2*DPTBV-AVDPT 

0255  IF(DAT( 10) ,GE .DMN2 . AND . DAT (10) .LE . DMX2)  DFRST=DAT( 10) 

0256  IF (DAT ( 10) . GE.DMN2 . AND . DAT (10) .LE . DMX2)  GO  TO  252 

0257  DLAST=DAT ( 10) 

0258  DPTBV-DLAST+AVDPT 

0259  GO  TO  20 

0260  C  SEE  IF  POINT  WITHIN  BOUNDS,  IFSO  ADD  TO  AVERAGE 

0261  33  DMX1-DPTBV 

0262  DMN 1 »DPTBV-AVDPT 

0263  IF (DAT( 10) ,GT. DMX1 .OR. DATC 10) .LT. DMN1 )G0  TO  250 
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0264 

0265 

0266 

0267 

0260 

0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

0277 

0270 

0279 

0200 

0201 

0202 

0283 

0284 

0285 

0286 

0297 

0288 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

0299 

0300 

0301 

0302 

0303 

0304 

0305 

0306 

0307 

0308 

0309 

0310 

0311 

0312 

0313 

0314 

0315 

0316 

0317 

0318 

0319 

0320 

0321 

0322 


IF (NUM8V. GT. 500) GO  TO  250 

C  DEPTHS  TOO  CLOSE  TOGETHER  CAUSES  ARITH.  PROBS  IN 

C  DIFFERENCE  METHOD,  SO  DON'T  ALLOW 

IF (DAT( 10) -DLAST.LT. 0.05) GO  TO  20 
252  XBV(NUMBV) =DENS 
DLAST-DATU0) 

NUMBV=NUMBV+1 
GO  TO  20 

C  AVERAGE  DENSITIES,  COMPUTE  QV,  AND  REINITIALIZE 

250  SUMRH--0 . 0 

C  DO  NOT  USE  INTERVAL  UNLESS  HAS  CONTAINED  AT  LEAST  4  VALID  PTS 

IF (IIUMBV.GT.5)  GO  TO  255 
NUM8V= 1 

DPTBV=DLAST+AVDPT 
GO  TO  20 

255  DO  260  J=1,NUMBV-1 
SUMRH-SUMP.H+XQVC  J) 

260  CONTINUE 

SUMRH=SUMRH/FLOAT(NUMBV- 1 ) 

CALL  BRVA 1 (DPTBV, SUMRH, RHO0,  BV) 

NUMBV= 1 

DPTBV=DLAST+AVDPT 
DAT (11) =BV 

C  FOR  PLOTTING  PURPOSES  PUT  DAT(10)  =■  AVG  FOR  INTERVAL 

DAT ( 10) = (DLAST+DFRST) /2 . 0 
WRITE (6, 270)DAT(10), SUMRH. DAT(II) 

270  FORMAT! IX, "DEPTH=".F10.2,4X, "RHO=",FI0.6,4X, "BV  FREQ  =  \F10.2) 

C 

C  TIME  PLOTS 

C 

756  IF(IAX(1) .NE. 12. AND. IAX(3) .NE. 12. AND. IAX(2) .NE. 12)  GO  TO  B4 
IF (NDAT1 . GT. I0UF2 (2) )GO  TO  15 

IF (NDHTl . EQ . IBUF2 (21 , AND .NDAT2 .GT. IBUF2(3) )GO  TO  15 

NDEL 1 “ IBUF2 (2) -HDAT1 

NDEL2* IBUF2 (3) -NDAT2+(NDEL 1*3600) 

DAT ( 12) * (FLOAT (NDEL2) /60. )  +  ( (FLOAT ( I-B) )/5760. ) 

IF(IAXd)  .EQ.  12) TX=X1AXL*X1 1  NCR 
I F ( IAX(3) .EQ. 12)TX=X2AXL*X2INCR 
IF ( IAX(2) .EQ. 12) TX“YAXL*YINCR 
IF(DAT ( 12) . GT. TX)  GO  TO  99 
C 

C  COMPUTE  DEL (T) /DEL (P)  AND  DEL (S) /DEL (P) 

C 

84  IF( !AX( 1) .EQ. 13. OR. I AX (2) .EQ. 13. OR. IAX(3) .EQ. 13) GO  TO  85 
IF( IAX( 1) .EQ. 14. OR. IAX(2) .EQ. 14. OR. IAX(3) .EQ. 14)GO  TO  85 
GO  TO  23 

85  DTMP-DAT(3)-TOLD 
DSAL  =D AT  d ) -SOLD 
DPRS=DAT (2) -POLD 

C  WRITE (6, 977) DAT(3) , DAT(l) . DAT(2) , DTMP, DSAL, DPRS 
C977  FORMAT ( IX, 6 (F 12 . 4, 2X) ) 

IF ( ABS(DPRS) . LT. 0. 00 1 )GO  TO  20 
DAT( 13) =DTMP/DPRS 
DAT (14)  -“DSAL/DPRS 
C  WRITE (6, 978) DAT (13)  ,DAT (14) 

C978  FORMAT ( 1X,2(F12.4,2X) ) 

T0LD=DAT(3) 

SOLD=DAT( 1) 

IF(POLD.NE.0.0) ITMP5=0 
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0323 

If (PGLD.EQ.0.0) 1TMP5-1 

0324 

P0LD*DAT  <2) 

0325 

IP ( ITMP5.EQ.  DGO  TO  20 

0326 

C 

032? 

C 

PLOT  OUT  THE  DATA 

0328 

C 

0329 

23 

CALL  PLTCDATC  IAXC 1) ) . DATC IAXC2) ) . DAT ( IAXC3) ) ,  IOPT) 

0330 

20 

CONTINUE 

0331 

GO  TO  15 

0332 

99 

ICOPY= 15427B 

0333 

WRITE (7, 754) ICOPY 

0334 

754 

FORMAT (A2) 

0335 

END 

0336 

SUBROUTINE  PLT(X1,Y,X2,  IOPT) 

033? 

COMMON  IREC, ISTOP, ISCAN,NSCAN, IFLOR 

0338 

COMMON  1BUF 1 (456) , IBUF2C456) 

0339 

COMMON  C1,C2,C3,C4, C5,C6, INSTR 

0340 

COMMON  ISHIPC6), 1ST, IYR, IMO, IDAY.LOU 

0341 

COMMON  X1AXL, XI IV, XI INCR, YAXL, YIV, YINCR,X2AXL,X2 IV.  X2 1 NCR 

0342 

X1MAX=X1 IV+CX1 INCR*X1AXL) 

0343 

YEXT=YIV+(YINCR*YAXL) 

0344 

YMAX=AMAX1 ( YEXT, YI V) 

0345 

YMIN=AMIN 1 (YEXT, YIV) 

0346 

X2MAX=X2 I V+ ( X2 I NCR*X2 AXL ) 

0347 

IFCX1.LT. XI IV. OR. XI .GT.X1MAX) RETURN 

0340 

IFCI0PT.EQ.2)  GO  TO  118 

0349 

IF (X2 . LT. X2 I V. OR . X2 . GT . X2MAX) RETURN 

0350 

118 

IFCY.LT.YMIN. OR. Y.GT.YMAX) RETURN 

0351 

XISCAL-XIAXL/CXIMAX-Xl IV) 

0352 

X1PLT=XISCAL*CX1-X1IV) 

0353 

X2SCAL=X2AXLAX2MAX-X2IV) 

0354 

X2PLT=X2SCAL*(X2-X2 IV) 

0355 

YSC AL = YAXL^ ( YMAX-YM I N ) 

0356 

IF  ( YINCR .  GE  .0)  YPl..T=YSCAL*(YMAX-Y) 

035? 

IF ( YINCR . LT. 0) YPLT=YSCAL*( Y-YMIN) 

0358 

NUMB  1=3 

0359 

NUMB2= 1 

0360 

CALL  SYMB (XIPLT,  YPLT, .  04,  NUI1B 1 » 0 . 0,  - 1 ) 

0361 

IFC IOPT. Ed. 2)  GO  TO  117 

0362 

CALL  SYMB (X2PLT, YPLT, .04,NUMB2,0.0,-1) 

0363 

11? 

RETURN 

0364 

END 

0365 

SUBROUTINE  INIT1 

0366 

COMMON  IREC, ISTOP, ISCAN,NSCAN, IFLOR 

036? 

COMMON  IBUF1C456), IBUF2C456) 

0368 

COMMON  Cl, C2,C3,C4,C5,C6, INSTR 

0369 

COMMON  ISHIPC6), 1ST, IYR, IMO, I DAY, LOU 

0370 

COMMON  X 1 AXL , X 1 I V, X 1 I NCR , YAXL , Y I V, Y I NCR , X2 AXL , X2 I V, X2 I NCR 

0371 

INSTR= 1 

0372 

TCTP=0.0 

0373 

C  1=0 . 1 

0374 

C2=O.0 

0375 

C3=O.0 

0376 

C4= . 0005 

0377 

C5=0 . 00 10 

0378 

C6=0.0 

0379 

IHTCP1 =4 

0380 

LAG=0 

0381 

URITEC1, 18) 

0382 

18 

FORMAT ( “USE  DEFAULT  VALUES,  1-YES  2-N0") 

0383 

READC l,*)KANS 

0384 

IF (KANS.EG. 1 ) GO  TO  99 

0385 

UR  I TEC  1,20) 

0386 

20 

FORMAT! "DEFAULT  VALUES  GIVEN  IN  PAREHTHS") 

0387 

UR  I TEC  1,22) 

0388 

22 

FORMAT C “ INSTRUMENT  *  Cl>") 

0389 

READ! I,*) INSTR 

0390 

URITEC 1,24) 

0391 

24 

FORMAT C "TEMP  PROBE  TIME  C0NST,SEC0NDS  (0.0)") 

0392 

READC 1,*)TCTP 

0393 

URITEC 1,26) 

0394 

26 

FORMAT! "CAL  CONST  PRESSURE  C1,C2  (0. 1,0.0)") 

0395 

READC  l,HOCl,C2 

0396 

URITEC 1,28) 

0397 

28 

FORMAT! "CAL  CONST  TEMP  C3,C4  (0.0,0.0005)") 

0398 

READC 1,*)C3,C4 

0399 

URITEC 1,30) 

0400 

30 

FORMAT! "CAL  CONST  COND  C5,C6  (0.001,0.0)") 

0401 

READC 1,*)C5,C6 

0402 

URITEC 1,40) 

0403 

40 

FORMAT ( “LAG  OF  TEMP  PROBE  #0F  SCANS  (0)") 

0404 

READ (1.*) LAG 

0405 

C 

UR  IT  :  OUT  VALUES 

0406 

URITEC 1,64) INSTR,  TCTP 

0407 

64 

FORMAT: IX. "INSTR  12, "TCTP  ",F8.4) 

0408 

URITEC i,66)Cl,C2,C3,C4,  C5.C6 

0409 

66 

FORMAT 'IX, “CAL  CONSTS  “.6F8.4) 

0410 

URITEC 1,68)X1AXL, XI IV.Xl INCR, VAXL, YIV,YINCR,X20XL.X2IV,X2INCR 

0411 

68 

FORMAT!  IX,  "AXES  CONSTS  \9F6.2) 

0412 

URITEC 1,70) INTCP1.LAG 

0413 

70 

FORMAT! IX, "PLT  TM  INT  ",I4, "LAG  IN  SCANS  ",I4) 

0414 

99 

RETURN 

0415 

END 

0416 

ENDS 
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FTN4.L 

SUBROUTINE  LRSAL(G.P.T.S) 

C  CTD  PACKAGE  SUBROUTINE  TO  OBTAIN  SALINITY  FROM  C.T.P 

C  MODIFIED  TO  FIT  INTO  MY  PLOT  PACKAGE.  18JULY77  L.  ROSENBLUM 

COMMON  IREC, ISTOP, ISCAN.NSCAN, IFLOR 
COMMON  IBUF1C456), IBUF2(456) 

COMMON  C1,C2,C3,C4,C5,C6, INSTR 
COMMON  ISHIP(6). 1ST, IYR, IMO. IDAY.LOU 
COMMON  PAXL.PIV. PINCR.TAXL.TIV, TINCR.SAXL.SIV, SINCR 
C  DEG  CELSIUS;  PRESSURE  IN  DECIBARS 
IF(ABS(S-35.0)-35.0)  3,3.2 

2  S-35.0 

C  SCALE  CTD  TEMPERATURE  TO  THE  1948  STANDARD 

3  CONTINUE 

XT-T+4.4E-6*T*( 100. -T) 

XC-G 

XP-P 

XTO  -  XT- 15. 

RC-XC/42.909 

AL  -  -7 . 6E-5*XT+6 . 95E-3 

AJ  »  <(-1.657E-4*XT+8.276E-3)*XT-1.535E-l)*XT+1.0 
RT  •  (<(-.67249142E-8*XT-.19426015E-6>*XT+.998B65B5E-4)*XT 
1+. 2013 1661E-1)*XT+. 67652453 

GO”  ( ( -7 . 9E-6*XT+8 . 3089E-4) *XT-4 . 5302E-2) *XT+1 .5192 
H- (-2 . 492E-9*XP+2 . 577E-5) *XP+4. 0E-4 
F«((3.3E-13*XP-3.3913E-B)*XP+1.042E-3)*XP 
SP  «GA*F  +H*AJ 
RCT-RC/RT 

RP0  ■  1.  +1 .0E-2#SP 
RP1-1.0E-2*SP*AL 
5  SO-S 

RP-RP0+RP1*(35.0-S) 

RS-RCT7RP 

R15  •  ((-.21E-5*RS*RS-.63E-5)*XT0+C37.3E-5*RS-72.E-5)*RS 
1+96.7E-5)*XTQ*RS#(RS-1  .0)  +  RS 

S-( (((-1. 3231 1*R 15+5. 98624)*R15- 10. 67869) *R 15+ 12. 80032) 

1*R 15+28. 2972) *R 15-0. 0B996 
IF (ABS(S-S0)-5E-2)  10,10,5 
10  CONTINUE 
END 
ENDS 
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FTN4,L 

SUBROUTINE  POTMPCPRS, TEMP, SAL, PTMP) 

C  POTENTIAL  TEMPERATURE  ACCORDING  TO  BRYDEN 
C  CONVERT  TO  1940  TEMPERATURE  SCALE 

TMP  *  TEMP  +  4. 4E-6*TEfP*( 100. -TEMP) 
S0-SAL-35. 

A3* .  50484E- 14*TMP-0  .  16056E-12 
A2» ( , 2 1 987E- 1 l*TMP-0 . 3 1 628E-9) *TMP+. 89309E-8 
1-.41057E-10*S0 
All*-. 29778E-7*TMP+. 17439E-5 
A 10- <  < . 40274E-9*TMP- . 54065E-7)*TMP 
1+.8319BE-5) *TMP+. 36504E-4 
Al-Al 1*S0+A10 

PTMP*TEMP- ( ( A3*PRS+A2 ) *PRS+A 1 ) *PRS 

RETURN 

END 


FTN4.L 

SUBROUTINE  SVEL(P,T,S,  V) 

C  J  UEBSTER  SVEL  SUBROUTINE 

C  COMPUTES  VELOCITY  OF  SOUND  IN  METERS/SEC  FROM  PRESSURE  IN  DECIBARS 

C  TEMPERATURES  IN  DEGREES  C  AND  SALINITY  IN  PARTS  PER  THOUSAND 

C  CHANGE  PRESSURE  TO  TOTAL  PRESSURE  IN  KG/CM*^ 

PKC  *  P  *  .1019716  +  1.03323 
SM  -  S  -  35.0 

C00  *  ( ( (7.985 1E-6*T  -  2.6045E-4)*T  -  4.4532E-2)*T  +  4.5721  )*T 

X  +  1449.14 

C01  *  (7.771 IE-7*T  -  1.1244E-2)#T  +  1.39799 

C10  -  ((4.5283E-8*T  +  7.48I2E-6)*T  -  1.8607E-4)*T  +  1.60272E-1 

Cll  ■  (1.5790E-9*T  +  3. 1580E-8)*T  +  7.7016E-5 

C20  -  (1.0563E-9*T  -  2.5294E-7)*T  +  I.0268E-5 

B3  -  -l.9646E-10*T  +  3.52I6E-9 

B0  -  (1.69202E-3*SM  +  C0l)*SM  +  C00 

B1  -  Cl 1#SM  +  C10 

B2  -  - 1 . 2943E-?#SM  +  C20 

V  -  (  (  (-3 .3603E- 12  *  PKC  +  B3)>KPKC+B2)*PKC+B1)*PKC+B0 

RETURN 

END 

END* 
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FTN4.L 

C 

SUBROUTINE  ATG(PRS, TEMP, SAL. GAMMA) 

C  COMPUTES  ADIABATIC  TEMPERATURE  GRADIENT  IN  DEG  CELSIUS/1000  DB 
C  CONVERT  TO  1948  TEMPERATURE  SCALE 

TMP  -  TEMP  +  4. 4E-6*TEMP*( 100.-TEPP) 

S  ■  SAL  -  35. 

A01  ■  -.42393E-4*TMP+. 18932E-2 

A00  •  ((.66228E-6*TMP-.6836E-4)*TMP+.85258E-2)*TMP+.35B03E-1 

A10  -  ( (- . 5448  IE- 10*TMP+. 8733E-8) *TMP- . 67795E-6) *TMP+. 1874 IE-4 

All  -  . 27759E-B*TMP~.  1 1351E-6 

A2  -  (-.216B7E-12*TMP+. 18676E-10)*TMP-.46206E-9 

A1  -  A11*S  +  A10 

A0  -  A01*S  +  A00 

GAMMA  -  (A2*PRS+A l) *PRS  +  A0 

RETURN 

END 

END* 


FTN4.L 

C 

SUBROUTINE  SIGMT(TrF.XS.XS0.XST) 

DOUBLE  PRECISION  XT.XS0.E1 .B1.B2.E2.XST 
C  CONVERT  TO  1948  TEMPERATURE  SCALE 

XT  -  TMP  +  4.4D-6*TMP*(100.-TMP) 

XS0  -  ((6.767861 36D-6#XS-4 . 82496 1 40D-4) *XS+8 . 1 4876577D- 1 ) *XS 
X  -9 . 34458632D-2 

El  -  ( ( (- 1 .4380306 lD-7n<XT-l .98248399D-3)*XT-5. 459391 1 1D-I)*XT 
X  +4.531 6B426D0 ) #XT 

B 1  -  ((-1.0843D-6n<XT+9.8185D-5)>t<XT-4.7B67D-3)*XT  +  1.0 

B2  -  ((1.667D-8*XT-8. 164D-7)*XT  +  1.803D-5)*XT 

E2-(B2  #XS0+B1)  *  XS0 

XST  -  E1/(XT+6.726D1  )+E2 

RETURN 

END 

END* 
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FTN4,L 

C 

SUBROUTINE  SPVOL(PRS,TErP.XSO,XST,XALPH) 

DOUBLE  PRECISION  A0,E2,XALPH.XSO.XST,TMP.BIO,B1 1.B12,B20,821 
1  B22.C3.C1.C2.E1 

C 

C  SPECIFIC  VOLUME  IN  CM**3/GM 
C 

C  CONVERT  TO  1948  TEMPERATURE  SCALE 

TMP  ■  TEMP  +  4. 4E-6*TEMP*( 100 . -TEMP) 

BIO-((4.0E-l2*TMP-6.63E-10)*TMP+3.673E-8)*T7P-2.2072E-7 
B 1 1 • (4. 0E- 12*TMP-3 . 28E- 1 0) *TMP+ 1 . 725E-8 
812-1. 0E-12#TMP-4.5E-11 

B20-C2. 14E-14*TMP-1 .  24064E-12)*TT'P-6.68E-14 

82 1 -  (-2 . 0E- 16*TMP+1 . 206E- 14) *TMP-4. 248E- 13 

822—  6.0E-17*TMP+t.8E-l5 
C3-1 .5E-17*TMP 
C1-(B12*XS0+61 1)*XS0+8I0 
C2- (B22*XSQ+B2 1 ) *XSO+B20 

E 1 - ( (C3*PRS+C2) *PRS+C 1 ) *PRS 

E2- 1 . 0-4. 886E-6*PRS/( 1.0+1. B3E-5*PRS) +€ 1 

A0- 1.0/(1 .8+.001*XST) 

XALPH*A0»E2 

RETURN 

END 

END* 


FTN4.L 

SUBROUTINE  BRVA 1 (DPT, RHO, RHO0, BV) 

DOUBLE  PRECISION  RDIF.E,RHO,RHOL,DDIF»GSPVO 
C  SUBROUTINE  COMPUTES  BRUNT-VAISALA  FREO.  IN  CYCLES/HR 

C  USING  APPROX.  GIVEN  IN  PHILLIPS  P.  17 

C 

C  PROGRAMMER  L.  ROSENBLUM  10/78 

C 

C  CONST  IS  GRAVITY  IN  M/SEC/SEC 

GSPVO-9.81/RHO0 
DDIF-DPT-DPTL 
RDIF-RHO-RHOL 
E«GSPVO*RD IF/DD IF 

C  CONST  CONVERTS  FROM  RADIANS/SEC  TO  CYCLES/HR 

BW-572 . 9578*DSQRT(DABS(E) ) 

BV-SIGN(BW.E) 

RHOL-RHO 

DPTL-DPT 

RETURN 


